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Multiple-Beam Interferometry : 
Intensity Distribution in the Reflected System 


Bye]. “HOLDEN 
Royal Holloway College, University of London, Englefield Green, Surrey 


Communicated by S. Tolansky; MS. received 29th October 1948 


ABSTRACT. The intensity distribution within the reflected fringe system from a 
multiple-beam interferometer employing silvered reflecting surfaces is discussed in its 
dependence upon the reflectivities and phase conditions at the reflecting surfaces. Using 
the Fizeau fringes of equal thickness localized on the zero order Feussner surface, the 
reflected intensity distribution is examined experimentally as the reflection coefficients of the 
interferometer surfaces vary between 4% and 90%. In the range of low reflectivities, 
when the silverings used showed marked colours, the changes in the reflected system lead 
to the measurement of a phase quantity related to the optical constants of the silver in the 
form of the thin film. Also in this range, the reflected fringes have a symmetrical form which 
is of use in the examination of sources of low intensity. In the range of high reflectivities, 
the conditions for the use of the reflected fringes in topographical investigation are discussed.. 
The findings apply to fringes of equal chromatic order as well as to Fizeau fringes in reflection. 


§1. INTRODUCTION 


HE reflected interference fringe systems find many applications in metrology, 
the “Fizeau”’ localized fringes of equal thickness being most frequently 
used. The surfaces of the interferometer usually carry a thin, partially 

reflecting layer of silver, but a general description of the reflected intensity distri- 
bution from such reflectors for the particular reflecting coefficients involved. 
appears to be confined to a little-known early paper by Hamy (1906), who used 
chemically deposited silver to form the reflecting surfaces. 

In the work to be described here, an investigation is made of the reflected 
intensity distribution given by an interferometer whose reflecting surfaces are 
formed by the thermal deposition of silver upon optical flats, so as to cover the 
range of reflectivities from 4% to 92°%. Experimentally, the Fizeau fringes 
localized on the zero order Feussner surface are used for the observations, but 
the results obtained also apply to fringes of equal inclination and to fringes of 
equal chromatic order (‘Tolansky 1945) in reflection. 

The observations of Hamy have been confirmed and extended. In the range 
of low reflectivities, when the reflecting silver films are less than 1004. thick, 
the relation of the measurable phase quantities associated with the reflected 
fringes to the more usual measurements of phase quantities in single films is 
discussed. Inthe range of high reflectivities the conditions for the use of reflected 
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fringes in topographical measurements, such as have been carried out using the 
transmitted system in the case of a diamond surface (Tolansky and Wilcock 1947), 
are derived. 


§2. THE INTENSITY REFLECTED FROM A PLANE PARALLEL SHEET 
The case of a plane parallel air film between silvered optical flats will first 
be considered. he properties of such a system as developed theoretically can 
be tested experimentally using a simple wedge interferometer if the wedge angle 
and separation are sufficiently small (Brossel 1947). 
Figure 1 illustrates two such reflecting planes, with normally incident, 


Tx 


Rx 
Figure 1. 


monochromatic, parallel light giving rise to families of reflected and transmitted 
wave fronts. 

Let the factors by which the amplitude of a plane wave is caused to change 
on reflection at the air—silver interface of the reflectors 1 and 2 be r, and 7, and 
let B,, By be the accompanying phase changes at these reflections. Let r,, « 
be the equivalent quantities for the reflection at the glass—silver interface of 
reflector 1. 

Finally, let ¢; and ¢, be the amplitude transmission factors of the surfaces 1 
and 2 and 7,, 7, the phase changes at these transmissions. The changes of phase 
B, Bo, « are illustrated in a geometrical manner in the figure, that is, by assuming 
them equivalent to a change in air gap, of arbitrary sign. 

Assuming a plane wave of unit amplitude to be incident normally upon th 
interferometer, the state of interference at any point X in the focal plane of a 
collecting lens is given by the vector sum of the disturbances due to each reflected 
wave. For the reflected system the resulting vibration at X is given by: 

Vx=r, cos(9—a)+ =X trkr”—} cos {[8 — 27,] —[(n—1)B, +n(B.+8)]}. 
1 


n= 
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Assuming that the absorption is the same for each direction of transit through a 
film then the reflected intensity is 


ple ae + 2tirsr,, cos (6— «+ By + 27,) — 2ttr37yr,, cos (27, — By —« 
+773 —2ryrz cos (d+ B, + Be) 
In the same nomenclature the transmitted intensity is 
tit3 
Tyg = 4.53 om = OB 
a i +7773 — 2172 cos (5+ B, + Bp) alfite) 


‘writing 


1 
5 Oo =D, 
{ize TEE S 
where A=5+[,+,, the optical separation of the two surfaces. This term is 
familiar in its réle in the Fabry—Perot interferometer, which makes use of the rapid 
variation of @, about 6+8,+8,=A=2nz7 as 47,>1. 


Since (27,—$£,—«)=F depends only on the phase conditions at the first 
Surface, 1, the reflected system now becomes 


Igy =72 + Onl (t2re)? + 2t¢ror,, cos (A + F) — 2t?r27,7, cos F]. 


For the simplified case of a symmetrical interferometer where r?=73=R and 
t?=t3=T, the expressions become 


Tex =O,lT"], 
Ip, =R+0,fT?R+2TR cos(A+ F)—2TR? cos F]. 


It is seen that while the transmitted system is described by an expression consisting 
of the function ©, (the Airy “‘shape”’ function), multiplied by a magnitude term 
- involving only the transmission factors of the two surfaces, the reflected system 
differs since in addition to a superposed background term FR there are two 
itransmission-like terms, one of which involves cos F in its magnitude; there is 
also a term which is unique, since it is out of phase with ©, by an amount F. 

The condition leading to maxima and minima in the expression for the 
reflected intensity Jp, 1s 


sin Af{tireryrs +1or,, cos F —r3rirer, cos F} 
-+cos.A{7,7, sin F+72r2rr, sin F} = 273r,r, sin F, 
‘which can be written 
msinA+ncosA=p, 


where m, nand p are functions of 7,, 72, 7,,f,and F. ‘The solutions of this equation 
‘can be written 
A=2n7 +a; A=(2n+1)r—-8), 


where one of the series of roots corresponds to maxima and the other to minima, 
the quantities a and b being of the order of 7 and dependent upon the values of 
14) To, Tz, ty and F, that is, upon the reflection coefficients of the two surfaces 
of the interferometer, the transmission coefficient and the phase changes at the 
first (incident) surface. The two series of roots will not, in general, occur at 
equal intervals. 

The conditions that the roots be equally spaced is that in the equation 
‘1 px/0A =0 the independent term p is zero, i.e. 2r3r,r, sin F=0, which can occur 

27-2, 
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when sin F=0. Two cases arise, namely, (i) F=2nz, and (ii) F=(2n+1)z- 
In the first case, intensity maxima occur at A=2n7 and minima at A=(2n + 1)z- 
In the second case minima occur at A=2n7 and maxima at A=(2n+1)z. 

The other limiting condition of the intensity distribution which arises from 
the form of the equation 0Z,,/8A=0 occurs when the two roots are coincident, 
that is when (iii) m?+n?=p? and 

22 
Fon tle ee ee 
29 Wo 

In the absence of any known relations between the phase quantities «, 81, 7, 
for a particular silver film it is seen that the reflected fringe system from an 
interferometer in which the silver film is the incident surface may be expected 
to be asymmetrical. Observation of this asymmetry with relation to the sym- 
metrical transmitted fringe system from the same interferometer will allow the 
empirical determination of F for different film thicknesses. Of the three limiting 
conditions when the reflected fringe system is symmetrical, it is found experi- 
mentally for films of silver upon glass that cases (i) and (iii) represent limits to 
which actual conditions can tend, while case (ii) arises when the glass is unsilvered. 


§3. EXTENSION TO THE CASE OF A WEDGE 

It has been shown in the case of a silvered wedge interferometer (Brossel 1947) 
that the phase of the mth transmitted beam with respect to the direct is closely 
given by 

py = Nd — n 8rte*/3r, 

as compared with ¢,,=78 for strictly parallel surfaces, where « is the wedge angle 
and ¢ the wedge separation at the point considered. It is found experimentally 
that if the terms n®87te?/3A are small the transmitted fringes are asymmetrical, 
and as the terms increase the asymmetry resolves into secondary maxima close 
to the main maxima. ‘The reflected fringes from a wedge also show asymmetry 
due to this effect but, in order to obviate confusion with the asymmetry due to.-. 
the intrinsic phase characteristics of surface 1, it can be made negligible by - 
choosing an appropriate wedge angle and order of interference to make the phase 
lag terms due to the presence of the wedge angle negligibly small. The deductions. 
made in the case of a parallel sheet thus apply to a wedge chosen in this way. 


§4. EXPERIMENTAL 

The optical arrangement used consists of a source arm complete with colli- 
mating system which can be moved with respect to a fixed wedge interferometer 
and viewing system, so that the transmitted and reflected fringe systems localized 
on the Feussner surface of zero order can be observed in turn. 

Special attention must be given to the collimation and the angle of incidence 
of the light in each of the two positions. It is found experimentally that the 
reflected intensity distributions close to the Feussner surface of zero order F, 
have a markedly asymmetrical appearance with respect to the appearance of the 
distribution on the Feussner surface. Focusing further away from the Feussner 
surface causes the asymmetry quickly to attain the appearance of an adjacent 
maximum and minimum against the general background. Figure 2 shows the 
appearance of the reflected system (below) when the plane of focus is relatively 
steeply inclined to the Feussner surface and intersecting it. The accompanying 
transmitted system photographed on the same plate (above) by means of a split 
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shutter camera, under the same conditions, does not reveal its asymmetry at normal 
exposures. Over a small field it is clear that this asymmetry will be confusing 
unless it can be made negligibly small. The distribution close to the surface 
F, has been studied (Brossel 1947) in transmission, and it obeys the phase relation 
$y = nd —n*47Xe?/r, where ¢,, is the phase of the nth beam relative to the direct 
at the point X, and where X is the distance in the direction normal to the wedge 
from the common line of the wedge of angle «. This relation has the same general 
properties as that described when considering the effects of wedge angle; it is 
therefore important to have X and « such that the terms n247Xe?/A are negligibly 
small. If the incident light is collimated so that the angles of incidence are 
between the normal and a maximum of 8, the position of the Feussner surface 
[Flo for the light incident at @ is (Figure 3) X,=(t/e) sin@, where X, is the 


ee Plane of Focus 


Figure 3. 


distance normal to the interferometer at a thickness ¢t. If the plane of focus 
coincides with the Feussner surface [Fo], for the normally incident light, it will 
be distant an amount {(¢/e) sin 6} from the surface [Fo|,, which is the Feussner 
surface for the light incident at an angle 6. Therefore 6 must be chosen such 
that the terms n747Xj«?/A are small compared with 76. 

Another reason why the collimation should be especially rigorous is that the 
incident beam has to illuminate the interferometer from two successive positions 
such that the transmission surface [Fo], and the reflection surface [Fo], are 
coincident within the same toleration as that discussed above for the Feussner 
surfaces within the one system. By using a collimating lens of focal length 
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20 cm. and a source aperture 1 mm. in diameter, @ is reduced to about 15’. Using 
a wedge angle of less than 1’, and working with an order of interference of less 
than 25, the focusing and wedge angle terms can be made negligible. 

Finally, the lens used for the observations must be of appropriate aperture to 
collect sufficient of the multiply reflectedemergent beams. The lens chosen was of 
focallength50mm.,f/3. Figure 2 (see Plate) shows the effects of using an inadequate 
lens. Apart fromthe curvature of the field, the point of interest interferometrically 
is that there exists a variation of fringe shape along the length of the fringe because 
the aperture of the lens is too small for the wedge angle used. In the figure 
are shown the transmitted (above) and reflected (below) fringes from the wedge; 
these errors are much more obvious in the reflected system where the light intensity 
is high and the background level between the fringes not close to zero as in 
the transmitted system where excessive over-exposure is needed to reveal 
them. 

The surfaces making up the interferometer were silvered in the following 
manner. The incident surface, 1, was prepared by depositing six thin layers 
of evaporated silver upon a clean glass flat so that each deposit was slightly dis- 
placed with respect to the previous one, overlapping it over most of its area. 
By this means thirteen small adjacent areas of different reflectivities were obtained. 
The extent of the entire deposit was such that the record of the fringes from all 
the areas could fall on one photographic plate. 

Preliminary experiments had indicated that a region of much interest was in 
the range of reflectivities of the first surface below 40°4. Consequently the layers 
of silver evaporated in each successive process were very thin, equivalent thickness 
about 30 a., so as to give adjacent areas of reflectivities differing by a few per cent 
over the range 4-60°%. ‘The second surface was a single silvering; three different 
silverings of 28, 58 and 87% reflectivities respectively were made for use 
successively as the back surface. 

The glass on to which the silver was evaporated was originally photographic 
plate selected for its flatness of surface over an area about 1 cm. square, and cleaned 
by immersion in a mixture of chromic and sulphuric acids for a period of days 
before being washed with distilled water, dried and subjected to ionic bombard- 
ment at a pressure of 10-2 mm. Hg in the silver evaporation plant. The evapora- 
tion of the silver was carried out at a pressure of 5 x 10-4+mm. Hg. 

The aperture exposing the area to be silvered could only be moved to its 
successive positions for each silvering by exposing the evaporation plant to atmos- 
pheric pressure, and thus allowing impurities to affect the surface of the silvering 
between deposits. ‘The silverings thus prepared were used in the interferometer 
within a few hours of their preparation, and no inconsistency in the observations 
due to possible changes in the films in this period was apparent. The thin 
silverings with reflection coefficients between 4 and 20° exhibited a striking 
range of colour both in transmission and in reflection, particularly from the glass— 
silver reflection, if the accompanying air—glass reflection was avoided. The 
colour changes were fully appreciated when the films of different thickness were 
matched side by side as in the case of these compound silvered interferometer 
surfaces. On one typical test surface six areas of progressively increasing 
reflectivity fellin the range below 20% reflectivity, and the film colours in trans- 
mission were clear yellow, yellow-orange, orange, purple, dark purple, blue. 
The films of reflectivity greater than 20°/, showed the characteristic blue trans-- 
mission colour of silver films. 
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§5. OBSERVATIONS AND MEASUREMENTS OF THE REFLECTIVITIES 
AN DEE PEAS EO U AN vei 

Figure 4 shows the localized Fizeau fringes from a wedge interferometer 
whose incident surface was formed of six silverings to give areas of reflectivities 4, 
7, 8, 11, 15, 23 and 58°, and whose back surface was formed of two silverings of 
reflectivities 28 and 58°, the field of view being arranged so that the seven areas 
on the front surface covered each of the different back surface areas. In the 
Figure the reflected and transmitted systems are indicated by R and T respectively, 
with appropriate subscripts to denote the reflectivities of the front and back surfaces 
for the particular area, the front surface being given first. In this particular 
compound silvering the final evaporation was rather heavy, bringing the maximum 
reflectivity of the area upon which it was deposited to 58%,; the second sequence 
of reflectivities on the surface was therefore made in the reflectivity range 35 to 
58°, which does not exhibit any large changes in F such as occur on the range 4 to: 
25°%,, and the fringes from these areas have been omitted in Figure 4. 

The figure shows the variation in intensity distribution of the fringes as 
the reflection coefficient of the first surface increases from 4 to 58°%. With the 
aid of the transmitted fringes, which have gradually sharpening maxima as the 
reflectivities increase from left to right, and of the two vertical and one horizontal 
lines which are the images of fine scratches ruled on the silver surface, it is easy 
to see the relative shift of the maxima and minima in the reflected system. When 
the first surface is uncoated glass (r?=0-04, F=(2n+1)z) the reflected system 
shows minima coincident with the transmitted maxima. (The area Roy. 5. 1S 
of some incidental interest since then the reflectors have greatly disparate reflec- 
tivities and yet the fringes are easily visible.) As the reflectivity of the first surface 
increases, the fringes go through an asymmetrical sequence until at about 15%, 
the reflected system shows sharp maxima almost coincident with the transmission 
maxima but much more luminous (the relative times of exposure of the two 
systems was 5:1). These fringes were first observed by Hamy (1906) and, 
noting their sharpness and high light efficiency, he suggested their use in the 
examination of weak sources. He gave only a verbal description of the fringes 
and little attention seems to have been given to his experimental findings. ‘The 
sharpness attainable by the “transmission-like”’ fringes is shown in Figure 5 (a), 
area R,;.g,, wherein the back surface is made to have a reflectivity of almost 
90°, and the light losses. are due only to the double transit through the 15°% 
reflecting silver film. 

As the reflectivity on the first surface is further increased, the fringes again 
become asymmetrical, the background intensity grows, and finally the visual 
appearance is that of fine minima. Above 65° a microphotometer is needed 
to detect the asymmetry, and above 85%, this asymmetry due to the phase relations 
at the first surface can no longer be certainly differentiated from that of other 
sources such as the wedge angle effect and the focusing defects. 

Since in the illustrations the reflected fringes were in each case recorded on 
the same plate as were the transmitted fringes, the relative intensity levels within 
the systems from each area are roughly comparable in Figure 4. Figure 5 (a) 
shows fringes from an interferometer made up of the same front surface as the 
upper right plate in Figure 4, but with a back surface of 87% reflectivity. It 
is seen that the only effect of increasing the reflectivity of the back surface is to 
sharpen up all the fringes of the sequence while they retain their relative shapes, 
which are derived from the tront surface conditions. 
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The transmitted and reflected systems shown in Figure 4 were microphoto- 
metered using a photographic recording instrument. Using the known plate 
response characteristic, the true fringe profiles were obtained, from which values 
of A satisfying the equation @/,,/@A\=0 were measured. Substitution of these 
values in the equation gives values of F(=2r7,—8,—«). ‘These are plotted in 
Figure 6 against the reflection coefficients for the air-silver reflections at the 
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front surfaces. ‘These reflectivities were measured by a Strong type reflectometer 
and confirmed by a graphical method from the transmission fringe shape, obtained 
by using the microphotometer in conjunction with the equation expressing the 
intensity distribution in the transmitted system. 


§6. THE REFLECTED SYSTEM AT LOW REFLECTION 
COEFFICIENTS 


The points in Figure 6 show a jump of 27 in the value of F in the region of 
15 to 20% reflectivity of the silver film. In the many reported investigations 
of the phase quantities « and £ of thin silver films some observers (Rouard 1937) 
have noticed a jump in the values of « at very small thicknesses, but the experi- 
mental difficulties (poor sharpness and visibility of Fabry—Perot fringes in trans- 
mission at the low reflectivities of the thin films involved) are considerable. 
Observation of the reflected fringe system from an interferometer, the front 
(incident) surface of which is the silver film under examination, immediately 
reveals the variation of the phase quantity F. From the discussion of symmetry 
of the reflected fringes, this experimental observation, that the fringes change 
from symmetrical minima when the front surface is unsilvered through an 
asymmetrical minimum—maximum form until they are almost symmetrical 
maxima when the surface is silvered to about 15°, means that F has changed 
from (2n+1)7 to 2mm. As the reflectivity further increases, the fringes first 
have an asymmetrical maximum-—minimum form and then asymmetrical minimum 
as the reflectivity increases to 60°, the minimum becoming sharper as the very 
highest reflectivities, 90°%, are approached. ‘This implies that F is tending to 
(2n +1)z, approaching it, however, from values between (27 +1)z and 2(n+1)z 
as the reversed sense of the asymmetry requires. Figure 4 shows the sequence 
for the wavelength 5461 a. 

The rapid variation of F in the range of very low reflectivities, and the colours 
which the metal films present, appear to be related to the state of aggregation 
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of the metal of the film, since these phenomena occur for equivalent film thicknesses 
which are very small for explanations based on interference. The measurement 
of F can be extended to all visible wavelengths if the fringes of equal chromatic 
order (Tolansky 1945) are observed with the interferometer. 

The relatively sharp ‘‘transmission-like”’ reflection fringes (cf. Figure 5 (a), 
-area R,;.7), appear to be of value if sources of low intensity only are available, 
since the reflected light only passes through a 15% reflecting silvering twice, 
-and is thus not seriously reduced in intensity. These symmetrical fringes occur 
when F=2nz, which condition, in the case of the wavelength 5461 a. and air— 
silver-glass, arises when the reflectivity air-silver is about 15°% (equivalent 
film thickness approximately 50a.). In the case of other metal films and other 
dielectrics this condition, fF=2nz, should occur at characteristic values of the 
film thickness. ‘This may mean that a greater reflectivity is possible on the 
incident surface of an interferometer, the reflected fringes still appearing as 
symmetrical maxima on a dark background. In the case of the system air— 
-silver-glass the ‘‘transmission-like” fringes have a half-width which is j the 
distance between orders, when the back surface of the interferometer is opaque 
(for the green region). These fringes are also of use for the interferometric 
-study of vibrating metallic surfaces. 


$7. THE REFLECTED SYSTEM AT HIGH REFLECTION 
COEFFICIENTS 

When the reflection coefficient of the first surface exceeds 60%, the reflected 
fringes present the visual appearance of fine minima on a background (Figure 5 (4)), 
-and a microphotometer is needed to detect the slight asymmetry still present. 
The present measurement of F in this region is not capable of high accuracy since 
‘the quantities to be measured are then small (shifts of A,,,, and A,,;,) and in the 
equation for the condition 0/,,/0A =0; asr7—1, 73 1, small errors inthe measured 
-quantities cause large variations in the calculated values of F. In this range F 
-can only be measured to an accuracy of 4/30, and for 0:-6<R,<1-0 is found to 
be equal to (2n +1)7 +7/10. 

From the earlier discussion it was seen that if the two conditions 
F— (2n+1)z, R,—1 could occur together the maxima and minima in the reflected 
intensity distribution would tend to be superimposed. It is observed experi- 
mentally that the reflected fringes (sharp minima) at high reflection coefficients 
become sharper than the accompanying transmitted system, whose sharpness is 
-controlled by the function ©, alone. It is thus apparent that in the reflected 
system the superposition effect consequent upon F—(2n-+1)7 is taking place. 
For the present investigation in the region of high reflectivities the above accuracy 
of measurement of F is sufficient. 

In order to examine the sharpness and contrast of the reflection fringes in 
‘the range 0:6<R,<1-0, which is the working range in topographical studies, 
we can take with good approximation F=(2n+1)z. Therefore the reflected 
intensity from a symmetrical, working interferometer possessing identical 
high-reflection coatings can be written (F=7; 77=73=R; 1] =t}=T) 


Ve pees T?R+2TR*—2TR cosA 
Rx = +4 1+ R?—2R cos A 


“The maxima occur at A=(2n+1)z and the minima at A= 2zn7. 
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Minima. The minimum intensity in the reflected pattern [Jp] min 18, therefore, 
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If there is no absorption in the silver films we can write R + T=1, andthe minimum 
reflected intensity in the case of no absorption, [Zgx]),in, becomes zero. If there 
is absorption defined by A, such that R+7+A=1, the minimum reflected 
intensity [Ipx]4i, becomes A?R/(1—R)?. Thus 


R : 
(eee = Velie (i —RYP A’. 


Therefore the minimum intensity in the reflected pattern is raised from its value 
zero in the case of no absorption to the value A?R/(1 — R)? when there is absorption 
defined by R+ T+A=1. 

Maxima. he maximum intensity in the reflected pattern is 


Lexlmax = + (do 2B ZY. 
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If there is no absorption we can write R+T=1, and the maximum reflected 
intensity becomes 
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If there is absorption defined as before, the .maximum reflected intensity for 
absorption A, [Tpx]*ay> 18 
R 


; roe 
Uaxlinax = R + (1+ R)? 


Since typical values of A at the higher reflectivities employed in multiple beam 
interferometry are of the order of ae A®<4A. Therefore 


[J xl peak =n 


{(3 -2R — R®) +(4?—4A)}. 
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Therefore when absorption A is present, defined by R+ 7 +A =1, the maximum 
reflected intensity is decreased by an amount 44R/(1 + R)?. 

The usefulness of transmitted fringes is usually described by reference to 
their “half-width” or the reciprocal of the half-width, which is a measure of 
sharpness, and to the quantity (J,.3,—Jmin}, which is a measure of their visibility. 
Similar indices can be used to describe the reflected system at high reflection 
coefhcients. The following Table (Tolansky 1946) gives typical values of 


k, T and A for silver films of quite low absorption as used in topographical studies 
(green light). 


FC) 60 70 US, 80 85 90 94 
Tl) 3525 27 fp) 16:5 10:5 a5) 0-7 
A(%) 4-5 3 3 35) 4°5 Bo) 593} 


Figure 7 shows the values of [Zpx|min and [Zpxlmax for these good practical 
values of R, T and A in the curves marked 1, while Figures 8 and 9 show respectively 
the ‘‘sharpness”’ and ‘“‘contrast”’ of the beneer and transmitted fringe systems 
from a symmetrical interferometer using these silverings. 
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i For comparison with reflected fringes it may be noted that the transmitted 
fringes from silverings of reflectivity 93 °/ have a half-width of 1 /40 buta maximum. 
intensity of only 5% of the incident intensity. 

Absorption. From the sharpness and contrast curves of the transmitted and 
reflected systems the superiority of the reflected fringes for precision measure- 
ments is apparent. The preparation of the silver films for use with these systems, , 


0 


100 


4100 


iG 
480 — 
i= 
= 
I es 
@ = 
= 405 NY 
(ey 
440 = B 
P=) 
= 
fon} 
20 , a 
! ! 0 : Jo 
60 70 80 90 100 60 70 80 90 100 
RY Rf 
Figure 8. Figure 9. 


however, requires extra care, as the following considerations will show. In 
Figure 7 are drawn two families of curves: the upper set shows the value of the 
maximum reflected intensity [Jpx|max for different values of the absorption A; 
the lower set shows the value of the minimum reflected intensity [Ix] min for the 
same values of the absorption. In each family the curve marked 0 shows the 
values of the quantities when A is very small, ~ 0-01, the curve marked 1 shows 
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‘the values of the quantities when A has the values shown in the Table above, 
.and the curves marked 2 and 3 show the quantities when the absorption has values 
respectively two and three times these good practical values. In the region 
of high coefficients (marked as H), it can be seen that a slight increase in the value 
of A for a film of reflectivity 90°% changes conditions critically from sharp fringes 
‘of good contrast to conditions in which the contrast is so poor that the fringes 
‘cannot be identified at all. For a film of an even higher value of R the increase 
in A required for effective fringe obliteration is quite small. This effect occurs 
in practice, and a silvering which immediately after preparation gives clear, sharp 
fringes can at times deteriorate in as little as 60 minutes to a condition giving 
markedly inferior fringes. It is as yet not possible to determine by eye whether 
a silvering will give a good performance in reflection, and it is therefore necessary 
to measure R, T and A with a simple reflectometer in order to calculate [Jpx|ihax 
and Velen: 

In the region of somewhat lower coefficients (marked as M), the effect of 
increasing the absorption is to decrease the maximum reflected intensity, while 
the corresponding increase of the minimum reflected intensity is negligible. 
This is illustrated by the fringes in Figure 5(c), where areas which have the same 
R but different A are shown. This was accomplished by cleaning a piece of 
glass chemically in a normal manner, then silvering part of its surface and stripping 
the silver off. The glass was again silvered, so as to cover the two differently 
cleaned areas to about 75° and then used as the front surface of an interferometer, 
the back surface being uniformly cleaned and silvered to 82%). The region of 
lower absorption is clearly seen in reflection, being characterized by the higher 
background intensity (hence better visibility); the minimum intensity was little 
changed in the fringes trom the two areas, as theory predicts. In the transmitted 
system the maximum intensity is less in the higher absorbing area, and the 
fringes are made to look sharper by relative under-exposure. The stripping 
process is thus seen to be an effective cleaning process. Areas of the same A 
but different R are recognized in Figure 5 (c) by the ‘‘step”’ between them. 

Collmation. When silverings of high reflectivity 90°, and low absorption 
are obtained (curves 1 in Figure 7), it is noteworthy that the reflected fringes are 
so sharp that the tolerances in collimation acceptable in the transmitted system 
must be decreased ; lack of attention to this means that the fringes are submerged 
in the background. 

Using the tolerances sigeeeens by Tolansky (1946), the following Table 
shows the source diameter D permissible (10cm. focal length collimating lens) 
when using the reflected system at R=90° for various interferometer separa- 
tions 7: 

t(mm.) 1 0-1 0-01 0-001 
D(mm.) 0:1 0°35 1:00 3-000 


§8. UNSILVERED BACK SURFACES 


Although it is desirable to have high reflecting coefficients within the limits 
‘discussed when using the reflected system in surface studies, cases arise when the 
back surface of the interferometer, which is the surface under examination, will 
not be coated with a highly reflecting film since the several techniques by which 
this is achieved may influence undesirably the surface structure. In this case, 
the natural reflectivity of the surface under examination will be used and will, 


Reflected Interference Fringes Aly 


in general, be considerably lower than that of silver. The effect of having R, 
small is to cause the shape factor @, to have less rapid variations about A=2znrz, 
and to decrease the value of the magnitude terms upon which it operates, in 
comparison with the background term R,, as the expression for the reflected. 
intensity from an asymmetrical interferometer (R,, R,) shows : 


es {R fs T?R,+2T,1/(R,R,) cos (A—F)—27T,R,R, cos “| 
a 3 1+R,R,—24/(R,R,) cos A 

The etfect of this is loss of contrast and sharpness. However, by decreasing R, 
to the range of values 60-70% the asymmetrical reflection fringes can be usefully 
employed since the sharp side of the fringe can be used for measurement while 
the contrast remains good. Even with R, as low as 25% the use of the- 
asymmetrical fringes will allow quite accurate measurements of topography to be 
made. Since the reflected fringes with R, greater than 60°% consist of some kind 
of minimum on a background it will be seen that any kind of effect which causes - 
the background to encroach on the fringe will give an apparent increased sharpness. 
For instance, in the range of asymmetrical fringes it is found by experiment 
and confirmed by theory that for typical silverings when R, is about 55% (a 
value dependent slightly upon the value of the back surface reflectivity R,), 
the minimum intensity within the fringes is close to zero. This allows a high 
degree of over-exposure when photographing the fringes and, although the fringes - 
are actually asymmetrical, the effect of the over-exposure is to obtain a final 
image resembling a simple minimum on a uniform background, of a sharpness 
somewhat resembling the fringes from 80% silverings. This is of use when 
dealing with back surfaces of low reflectivity. Other particular procedures . 
are of use when dealing with certain surfaces, and a report of the actual investiga- - 
tions on several metallic surfaces will be published later. 


§9. CONCLUSION 

From the point of view of topographical studies, the characteristics of the - 
reflected system may be summed up briefly. For reflectivities greater than 90°% | 
the fringes are sharper than the corresponding transmitted fringes but their 
visibility is critically dependent upon the absorption in the silver film; collimation 
conditions are strict. Complications of the simple interference relationship . 
due to focusing defect and large wedge angles are immediately obvious in the 
appearance of the fringes, and one can avoid the errors that occur in transmission 
work if high precision measurements be attempted without sufficient attention 
to these factors. Finally, the results given apply to the reflected fringes of 
equal chromatic order which obey the same fundamental formula. 
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ABSTRACT. ‘The change in optical transmission of a solution of colloidal graphite in 
a magnetic field has been measured as a function of a field strength and of concentration, 
and the results obtained compared with those to be expected from theory. 

It is found that the experimental results agree closely with the hypothesis that the 
change in transmission is due to the change in projected area of the plate-like particles 
when these are aligned by the magnetic field. 

By comparison of theoretical and experimental curves the particle volume is found. 
The rate of decay, after removal of the field, of the change in transmission produced by 
the field has also been measured. 

An approximate theory is given to connect this rate with the diameter of the particles, 
and from this the size of the particles has been calculated. The dimensions so found 
agree satisfactorily with evidence from other sources. 


§1. INTRODUCTION 

, HE change in optical transmission of a medium consisting of a suspension 

of small particles in a liquid produced by application of a magnetic field 

in the same direction as the light beam was first observed qualitatively 
by W. R. Groves in 1845 (McKeehan 1940). He used for his experiment a 
suspension of ferromagnetic particles in water. his effect has recently been 
‘the subject of several papers (Heaps 1940, Elmore 1941, Mueller and Shamos 
1942), but the theory remains incomplete owing to the doubt which exists about 
the magnetic properties of these very small ferromagnetic particles; nor has a 
-complete experimental investigation been carried out. 

A similar effect has been observed by Stuart (1938) in a solution of colloidal 
-graphite. In the course of experiments which he carried out to find a method 
of detecting and estimating small quantities of graphite in oil, Stuart found that 
the optical transmission of such a solution changed when a magnetic field was 
applied, and he made approximate measurements at a fixed field strength for 
various concentrations. 

As graphite has known magnetic properties, being diamagnetic and anisotropic, 
it is possible to develop a satisfactory theory of this effect. The object of this 
paper is to describe experimental results obtained, and to show what information 
about the size and shape of the graphite crystals can be obtained by comparing 
the theoretical and experimental results. 


2 ee O RY: 
(1) The Change in Optical Transmission as a Function of Field 
Strength and Concentration 
A solution of colloidal graphite in water forms a true colloid system, in that 
the particles remain in permanent suspension and exhibit Brownian motion. 
As the particles are disc-shaped, the absorption due to any one particle will 


Optical Transmission of a Colloidal Solution 419 


‘depend on its orientation with respect to the light beam, and may be assumed 
proportional to its projected area; the diffraction effects which will also enter 
into the total absorption may be neglected here, as they may be assumed inde- 
pendent of orientation, and we are concerned here only with that part of the absorp- 
tion which is dependent on orientation. 

In order to find the relation between absorption, or transmission, and the 
applied magnetic field strength, it is first necessary to find the distribution of 
orientation as a function of field strength. This is determined by two factors— 
the torque exerted on the particle by the magnetic field, and the disturbing effect 
of the Brownian motion, which acts to prevent perfect alignment, tending always 
to re-establish random orientation. 

Distribution of orientation as a function of field strength. ‘The orientation of a 
particle may be defined by two coordinates giving the position of the end of the 
normal to the basal plane on the surface of a unit sphere whose polar axis is 
parallel to the applied field. ‘The coordinates used here are: 6, co-latitude, 


e 


7=M Hsin@ 


SS 
Sf, = -X,Hc0s 8 


(M,= -X, Hsin @ 
73=M, cos 0 
Figure 1. 
and ¢, longitude. The torque acting on a particle of volume V whose sus- 
-ceptibilities are —x,, —xX2 normal to the basal plane and in the basal plane 
respectively is given by (see Figure 1) 
T=T,-T,=— x,VH cos H sin8+y,VH sin@ H cosé 
=(v2—x1)H?V sin 6 cos 8, 
where H is the magnetic field strength. The potential energy (rotational) is 
therefore given by 
U=fT d0=3(x, — x2)H?V cos? 8. 
From Boltzmann’s equation, the number of particles having orientation between 
6 and +86, and ¢ and $ +5¢ will therefore be 


oo (X41 exe) 2 2 a| : SO 5 
dN =C exp | - oT HV cos*@ | sin 6 60 64, 
which, after integrating to find C, gives | 
SN _ 2x exp (—2? cos? @) sin 6 36 d6 (1) 
FNCU RRIF MOM Soe esa cote: 


where =a (u— mV, and B= | e-" dey, 
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Light parallel to field. For a beam of light parallel to the applied field, the 
projected area of a particle of area a is a cos, assuming the thickness negligible 
compared with the diameter. ‘The total projected area of the N particles per 


unit volume is therefore 


nl2 m2 
ae | exp(—a? cos*6) cos sind dbd$. ...... (2) 
vA 0 /0 


Integrating this with respect to ¢ and 6, we have A, = Na{1 —exp(—«?)}/2«E,. 
For zero field strength, «—0, and, therefore, Ay) = Na/2, so that 


Ay= 


ee Ay PS exp (ar) 
ipo aa Ca 
Light normal to field. When the light beam passed at right angles to the applied . 


field, the projected area of a particle of area a is a sin 6 cos ¢. 
From (1), therefore, 


pe CNG tg pope ane . 
Ay= 7Eas |, exp (—«? cos?@) sin?@ cos¢d6 dd, ...... (4) 
which gives 
y 7/2 
mn = = | jeep (or cost) sintP db. | 2222: (5) 


Relation between my and change in optical transmission. ‘The optical trans- 
mission of a medium containing particles having total projected area A is easily 
shown, as in deriving Beer’s law, to be given by T=e~4. Hence, writing 7, for 
the transmission in zero field, and 7, for the transmission in a field H, we have 


che 
Fi =€xP(4o~ An) =exp (4g(1 —mz)} 


Hence In — =A,(1—mg)= a — My). 
the Z 


Summary. It has thus been shown that the change in transmission of a 
solution of colloidal graphite in a magnetic field is given by 


log (Pai ty) =KAS ig) ~ Se ee (6) 


where K is a constant, proportional to the concentration, and my is a function 
of the field strength, given for light parallel and normal to the applied field by 
equations (3) and (5) respectively. 

Graphs of (1 — my) against the parameter « are shown in Figure 2. 

Equation (6) also shows that for a constant field strength log(7¥4/T,) is 
proportional to concentration. 


(ii) Theory of Relaxation Time 


If the particles are aligned by a magnetic field, and this field is then removed 
suddenly, the Brownian motion of the particles will re-establish random orientation. 
The time required for the process will be a function of the particle radius. Thus 
the change in optical transmission produced by the field will decay until the. 
transmission has returned to its normal zero-field value. 

Debye (1929) has given a formula for the relaxation time for polar molecules; 
this formula is not, however, applicable to this case, for two reasons: first, the 
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dependence of the observed effect on angle of orientation is different, and second, 
Debye’s formula applies only to small degrees of alignment, i.e. cases in which 
the curve of distribution versus angle of orientation is broad. An approximate 
formula for high degrees of orientation is therefore derived here, following 
Einstein’s method for the case of diffusion from a point (Einstein 1906). 

At the moment of removing the field, the particles are aligned with their basal 
planes parallel to the field, and the coordinate @ defined above has the value 
37 for this position. Itis more convenient now to take a new coordinate p=tnr—-O. 


Log (1-7y) 


Field 1 to Light 


Figure 2. 


If the number of particles having orientation between ¢% and % + 8y at time f 
after removal of the field is 
ON = Fs, t) cospde, woo see Ly) 


then the change in the number of particles between and +6 in a further 
short time d¢ is 
OF 
o°N =8t vs COS Orrin, at a cence (8) 
If the number of particles suffering a change of orientation of value between A 
and A +6A in time 6¢ is 
dN = Nf(A)oA, 
then the number of particles entering the zone between ¢ and 4+8y% from the 
zone between +A and 4 +A+6A in the time 5¢ is given by 


oN = Fi’ +A, t) cos (ob + A)SAPS(A) Seb 
=cospF(y+A, z)f(A)6A approx.  ...... (9) 
cos (+A) being replaced by cos in the approximation, since, provided 5¢ is 
small, only small values of A need be considered. ‘This approximation is equi- 
valent to regrading each elementary zone of the spherical surface as a section of 


a cylinder. 
PROC. PHYS. SOC. LXII, 7—B 28 
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Expanding F( +A, ¢) and integrating over all values of A, we obtain for the 
total number of particles in the zone between ys and + dy at time ¢ + d¢ 


SN =cos yp 5yp | Fu, t) ee f(A)dA+ s = Af(A) dA 
12eF 


sae | Ons recite W| 


This automatically includes the particles in the zone , +6) at time ¢ which 
still remain there at time ¢+5t, as these correspond to values of A less than dx. 
Since 


emu] 2 on|2 
| f(A)dA=1, — and | Af(A) dA =0, 
—an/2 —n/2 


this reduces to 
Lar 2 
SN =cos pf Sys (FU, +555 | A2f(A) as), Pes: (10) 
pb —n/2 


so that the change in the number of particles in the zone between #% and J +d% 
in time 6¢ is 
1 oe Ree 
2 ae = 2 
BN = cos 5h 5 5 cy dA 


=cos hb 5 Sr OF 


Equating this to the value of 62 found in equation (8), we have 
= =sa5T RT BOO. Gane (11) 


since from Einstein’s equation for the rotational Brownian movement, 

3(A?/5t)=RTB, where R is Boltzmann’s constant, T the absolute temperature, 

and B the rotational mobility (i.e. angular velocity per unit torque). : 
The solution to equation (11) is 


Fis, t) = Ct-* exp (—?/B2), where B=4kRTB. 
Hence oN = Ct exp (—#?/Bt) cos pd 
BN __ exp (—¥2/t) cos dp 
n|2 2 
J, exp (—¥IBE) cos pa 


from which the projected area of the particles at time ¢ is given by 


7/2 
Na | sa A ae LIC CI cea | 


Ay = SS ae eae a a ee ee 
J exp(—#IBE) cos fay 2 


For large values of t, %,/%,—>4, so that A,>4Na, which is the normal value in 
zero field; calling this A,, we have m,=A,/A,.=24%,/%. The quantity 
m, is connected with the transmission by a similar relation to that given for my 
in §2 (1) above, i.e. log (T,/T..) = K(1—m,). 


and 
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In Figure 3, log(1—m,) is plotted as a function of Bt; from Bt=0-2 onward 
it approximates to a straight line, showing that (1—®m,) decreases exponentially, 
the mean time constant being given by Br=0-8, or r=0-8/8 approx. Hence 
it is to be expected that log(7;,/T.,) will decrease exponentially after the field is 
removed. Before this relation can be used to determine the particle radius R, 
it is necessary to fix the value of B, in which R is contained; for spherical particles 
B=1/87R*. Perrin (1934) has compared the values of B for ellipsoidal particles 
with those for spherical particles of the same volume. — For an ellipsoid of revolu- 
tion whose major and minor axes are in the ratio 5:1, he finds B about 2-5 greater 
than for a sphere of the same volume, hence about half the value for a sphere of 
radius equal to the semi-major axis. As we wish to use the relaxation time to 
find R, and R « B-?, an error of a factor of 2 in the value assumed for B will not 
affect the result by a factor of more than 2', i.e. 1:26; the error in using the value 
1/877R? will probably be less than this, as departures from Stokesian conditions 
will tend to make the actual value of B larger than the theory indicates. The 
walue 1/877? will therefore be used here. 


Log (1-/71¢) 


Figure 3. 


The time constant 7 will-then become 7=0:8 x 277R?/RT, which gives 
Re kRTr é 
1-677 
§3. EXPERIMENTAL METHODS AND RESULTS 
(i) Experimental Methods 

When the applied magnetic field is small, and the concentration low, the 
change in transmission may be only a small fraction of the transmission in zero 
field. It is necessary therefore to use a method of measuring the change in which 
the effect due to the magnetic field will not be masked by fluctuations in the 
intensity of the light source or in the sensitivity of the photocell. This is achieved 
by using a comparison method, in which the intensity of the light beam passed 
through a cell containing the specimen of colloidal graphite is compared with 


the intensity of another beam from the same source, the second beam being 
28-2 
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adjusted to be equal to that through the graphite in zero field. Two different 
methods of comparison have been used in these experiments; the first uses a 
single photocell to compare the two beams, which are switched on to it alternately 
by a rotating glass disc, half of which is silvered and half clear. If the two beams 
are unequal in intensity, an alternating voltage of square waveform appears across 
the photocell load resistance, of amplitude proportional to the difference of 
intensity between the two beams. This is amplified and applied to the Y plates 
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Figure 5. 


of a cathode-ray oscillograph. ‘The glass disc rotates at 50 r.p.s., and the oscillo- 
graph time-base is set to run at 100c/s., so that alternate half-cycles of the square 
wave appear on alternate traces, giving two parallel lines on the cathode-ray 
tube. When the intensity of the two beams is the same, the lines coincide. 
The intensity of the reference beam having been so adjusted with no magnetic 
field on the specimen, the change in transmission produced by switching on the 
magnetic field will produce the square-wave voltage, causing the two traces to 
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separate. ‘I'he amplitude of the voltage produced at the photocell is then 
measured by applying another square-wave voltage, in the opposite phase, to 
the lower end of the photocell load resistance, the amplitude of this voltage 
being adjusted to make the traces coincide again. This voltage is derived from a 
D.C. supply via a commutator driven on the same shaft as the glass disc, and it is 
the D.c. voltage which is actually measured. This method gives greater accuracy 
than can be obtained by attempting to measure the voltage of the original square 
wave directly, and also ensures that the voltage across the photocell is maintained 
constant, although the two beams are of different intensity. 


. Photocells 


Figure 6. 


Figure 7. 


The second method of comparison uses two photocells, one for each beam. 
‘The outputs from the two ceils are amplified by the two halves of a double triode, 
|. between the anodes of which is connected a milliammeter. The equality of the 
two beams is indicated by zero current on the milliammeter. The change in 

ntensity of the beam through the colloid specimen is measured by finding the 
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voltage which has to be applied to the lower end of one of the photocell load 
resistances to restore the balance. Thus in this method, as in the previous one,, 
the voltage across the photocells is kept constant. 

The optical and electrical arrangements of the first method are shown dia- 
grammatically in Figures 4 and 5, and of the second method in Figures 6 and 7. 


(ii) Experimental Results 


Of the two methods described above, the firstis rather more accurate, due to the 
absence of drift and slow fluctuations, which is the result of using an A.c. amplifier. 
This was therefore used to measure the change in optical transmission as a function 
of the applied field, for light parallel and normal to the field. The results obtained 
are shown graphically in Figure 8, where the experimental points are shown super- 
imposed on the theoretical curve, which has been drawn to such a scale that its 
saturation level coincides with the extrapolated saturation level of the experimental 


H (Oersteds) 
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Figure 8. 


points. The agreement is very close, apart from the length of the toe at low field 
strength. ‘This is less in the experimental curve than in the theoretical one, 
because of the distribution of particle sizes in the solution used; as the parameter 
« contains the product H?V it will readily be seen that large particles will apptoach 
saturation at much lower values of field strength than small ones; hence the 
effect of a small number of large particles present in the solution is to steepen the 
first part of the curve, thereby reducing the length of the toe. 

Effect of varying the concentration. 'The second experimental method described 
above was used to study the effect of varying the concentration. Although this 
method is not quite so accurate as the first one, it is more readily adapted for 
measuring the transmission factor 7; this was not required in the first experiment, 
as relative values of the change produced by the field were all that was required. 
In these experiments, however, accurate values of 7) are required since, the 
concentration being very small (of order 10-* by weight), more accurate results 
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may be obtained by measuring 7, and assuming the validity of Beer’s law than by 
measuring the concentration directly. 

The results are shown in Figures 9 and 10. In Figure 9, log(Ty/7T,) is plotted 
against H for each concentration; it will be seen that all the curves are similar, 
as is to be expected. In Figure 10, log(Ty/Ty) is plotted against log Ty 
(oc concentration) for various field strengths. There is some falling off from the 
direct proportionality predicted by the theory at higher concentrations; this 
may be due to increasing importance of diffraction effects. 
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Different particle sizes. A solution of colloidal graphite of concentration of 
order 10~-? by weight was allowed to sediment in a tall vessel, and then separated 
into anumber of layers. In Figure 11, log(7,/T ) is plotted against H for three 
of these layers, Layer 1 being the top, Layer 2 a central one, and Layer 3 the 
bottom. It is noticeable, particularly in the cases of Layers 1 and 2, that the 
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greater homogeneity of particle size in these specimens as compared with the 
normal solution has lengthened the toe of the curve, making it similar to that 
predicted by the theory. 

Calculation of particle volume from these results. The parameter « defined 
in §2(i) is [(v,—x»)/2RT}HV}, which is proportional to H. By comparing the 
experimental curves in Figures 8 and 11 with the theoretical curve for log (T/T ) 
we can obtain a value of the factor «/H=[(x1;—x2)/2RT}'V?. The values of k, 
Boltzmann’s constant, and 7’, the absolute temperature, are known, and Ganguli 
(1936) gives values for x, and y, for small crystals, of the order of size to be 
expected here. 

His values are — y,= —0:9 x 10-® per gm., — x,= —0:5 x 10-§ per gm., or 
X1—X2=4x10-7 per gm. The values of «/H, together with the values of particle 
volume calculated from them, are given in Table 1. 


Table 1 
: Volume 
Specimen o/H fom? x 10-22) 
Normal solution 1/450 0:44 
Layer 1 1/750 0:16 
Layer 2 1/600 0-25 
Layer 3 1/300 1-00 


(111) Relaxation Time Measurements 


Method of measurement. The second form of the apparatus described above 
was used to measure the relaxation time. The only change required was the 
removal of the milliammeter between the triode anodes; the voltage difference 
produced by the change in transmission due to the field was then amplified by a 
p.c. amplifier and applied to the Y plates of a cathode-ray oscillograph. The 
time-base was set to run at about 3-5 cycles per second, and on switching off 
the field the decay curve was seen spread out over several traces. This was 
photographed together with a timing trace, and from measurements taken from 
the photographs the four graphs in Figure 12 have been plotted. The quantity 
which is to be expected from the theory to decay exponentially is log(7,/T..), 
so the ordinate in Figure 12 is log log(7,/T.,), plotted against time in seconds. 

The time constants 7 found from Figure 12 for the normal solution and the 
three sedimented layers are respectively 0-59, 0:17, 0:29, 0:89 second. 

From these values, by means of the relation R =(R77/1-677)', we can calculate 
R for each of the four specimens. ‘Then from R and the values of V given in 
Table 1 the shape factor, the ratio of thickness to diameter of the particles, can 
be calculated. Table 2 gives the values of volume, diameter and shape factor 
for the four specimens. 


Table 2 
Specimen Mol Diameter (uw) Shape factor 
(cm? x 10-1?) s P 
Normal solution 0-44 1-56 1/6°9 
Layer 1 0-16 1-04 1/5-45 
Layer 2 0-25 1-24 1/5-95 


Layer 3 1-00 1-80 1/4-6 
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§4. A.C. MAGNETIC FIELD 

Since the torque acting on the particles is proportional to H?, it is to be expected 
that an alternating field will have a similar effect to a steady field. This was found 
to be the case, and in ‘Table 3 comparison is made between steady and alternating 
fields, the change in transmission here being expressed as a percentage. In 
the case of the alternating current, the current given is the peak value, not r.m.s 
It will be seen that at fairly large field strengths the change in transmission produced 
by the a.c. field is nearly 90°% of that produced by a steady field equal to the peak 
value of the alternating field. This is due to the fact that the relaxation time is 
much longer than the orientation time; the latter decreases with increasing field 
strength, whereas the former remains constant. 


‘Tables 
DSC! 

Current (amp.) Aompeak \ 0-25 0:5 1-0 ey 2:0 255 
(a) 5T/T, steady field 0:066 0-145 0:256 0:334 .0°381 0-413 
(b) 5T7/T, alternating field 0)-03207 0-038aun OF6Gs—8 022550 02325 nm 023.60 
(6) as percentage of (a) 49 61 65 76 84 87 


Note: 1 amp. corresponds to approximately 1,000 gauss. 


§5. EXAMINATION OF PARTICLES UNDER MICROSCOPE 


An independent check on the values found for the diameter and shape factor 
-of the particles in the experiments described above was obtained by examination 
under a microscope. The colloidal graphite was mixed with gelatine and allowed 
to set in a magnetic field, being slowly rotated about an axis normal to the particles 
with their basal planes all parallel; thus when the gelatine has set, sections can 
be cut normal to the basal planes. When such a section is viewed under an oil- 
immersion microscope, the particles are seen edge-on, presenting a rod-like shape. 
It is not of course possible to measure the shape factor of the smaller particles, 
of order 1 » or less in diameter, as these, seen edge-on, are at the limit of resolution ; 
approximate measurements on a few of the small proportion of larger particles 
(2-5 w) present gave values of about 1/5 for the ratio of thickness to diameter. 

The average particle diameter is known, from measurements carried out by 
Acheson Colloids Limited, to be about 1-1, which agrees sufficiently well with 
the results obtained in these experiments, when it is remembered that, in these 
optical absorption measurements, the larger particles present, although small 
in number, have a disproportionate influence on the mean value, owing to their 
large area. 


§6. CONCLUSION 


The results of the experiments carried out to find the change in optical trans- 
mission of the solution of colloidal graphite as a function of the applied magnetic 
field agree well with theory given in §1. This theory is based on the assumption 
that the change in optical transmission depends only on the change in projected 
area which occurs when the particles are aligned by the field. From the agreement 
found between theory and experiment it appears that this assumption is largely 
justified. 

The value found for the mean particle volume from these experiments 1s 
‘0:44 x 10-12cm?, and from measurement of the relaxation time the mean diameter 
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is found to be 1:56. Both these values may be expected to be larger than the 
averages by number found from direct microscopic measurement of the particles, 
as the larger particles, having greater area, will play a predominant part in deter- 
mining the change in optical transmission produced by alignment; thus the mean 
value will be displaced in the direction of larger particle size. "The value given 
for the average particle diameter by Acheson Colloids Limited, whose preparation 
Aquadag ‘S’ was used for these experiments, is 1-1. This agrees sufficiently 
well with the value 1:56 found here, when this displacement is taken into 
consideration. 

The value of about 7:1 found for the ratio of diameter to thickness appears 
reasonable, being of the same order as that usually found for macroscopic crystals 
of graphite. It is noticeable that this ratio remains approximately constant over 
the range of particle sizes examined here. 

From the results of these experiments it can be concluded that reasonably 
accurate information about the size and shape of colloid particles can be obtained 
by studying the change in optical transmission which takes place when the particles 
are aligned, at any rate in the case of colloidal graphite. It seems likely that 
similar information could be obtained for other colloids in which alignment 
effects can be observed, in magnetic or electric fields. 
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ABSTRACT. A general, synthetic method of obtaining rigorous solutions of steady 
electron flow subject to space-charge forces is presented. The solutions are not obtained 
for given boundary conditions, but the boundary conditions are deduced from the solutions. 
Two examples of such solutions, involving strongly curved two-dimensional electron 
trajectories, are given; the method is in principle capable of giving the solutions of all 
possible electron flow patterns in three dimensions except perhaps those involving inter- 
crossing trajectories. It is suggested that the subject offers scope for applied mathe- 
matical research at least on the same scale as potential theory. 


§1. INTRODUCTION 

s far as the writer is aware, no rigorous solutions for electron flow subject 
to space-charge forces have been published, except for cases involving 
straight line trajectories. The latter are involved in the classical solutions 
(Langmuir 1913, Langmuir and Blodgett 1924) for planar, cylindrical and spherical 
symmetry leading to the Child—Langmuir three-halves power laws, as well as in the 
more recent application of non-linear mechanics to the problem by Bellustin (1939). 
It is surprising that other types of motion do not appear to have been thoroughly 
studied, when one considers the complicated structures of modern electronic 
devices and how very far from rectilinear their electron trajectories often are. 
For example, the important question of current division between grid and anode 
in an oscillator triode can hardly be investigated seriously until the trajectories 
in the neighbourhood of the grid are known; only then can one take into account 

such factors as secondary emission from the grid. 

This paper presents a synthetic general method for obtaining rigorous 
solutions of steady electron flow. By “rigorous”? we mean that space-charge 
forces are taken into account; but, as in the original derivation of the Child— 
Langmuir laws, we shall neglect random (temperature) variations of velocity, 
assuming that the steady flow is characterized by the existence of a unique velocity 
vector at every point in the stream. Similarly, we shall assume that particle 
velocities are so much smaller than the velocity of light that forces due to the 
magnetic field of the currents may be neglected. 

After an exposition of the basic idea of the method, the results for two cases 
of curvilinear two-dimensional flow are presented, followed by a general discussion 
of the achievement of these and other solutions. 


§2. FUNDAMENTAL EQUATIONS AND METHOD OF SOLUTION 
The motion of an electron is governed by the force equation 
mv = —eE, were cL) 
where e, m are the magnitudes of the electronic charge and mass, and v, E are 
respectively the velocity of the electron and the electric field acting on it. 
If we consider the x component of the velocity in a right-handed Cartesian 
rectangular coordinate system with unit vectors I, J, K we have 


dvs) easieradaiee e188 Base. (2) 
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where ds is a differential displacement in the direction of the motion. Hence 
Us=VoTad Ue, == 5 « | | eenaete ot ae 
Similarly Dy =v Brad vy, lke) i ye sare (4) 
and C= Viotadv, |) = a). nen ee (5) 
Substituting these three expressions in equation (1), 
I(v gradv,) +J(v gradv,) +K(v gradv,)=—(e/m)E. ...... (6) 
Furthermore, the following relations have to be satisfied : 
Cty (oV) 0,0 deen oar neeenl senate (7) 
div ES pieje  * Ps ea es ere, (8) 
B= gradi) 6 fo tak ©) She Me Oe (9) 


where p, V are respectively the charge density and the potential and eq is the 
permittivity of free space (in rationalized M.K.s. units). These equations express 
respectively the continuity of current, Poisson’s law and the irrotational nature 
of the electric field. . 

Equations (6), (7), (8) and (9) are equivalent to eight simultaneous equations 
in eight unknowns, namely, three components of velocity, three components 
of electric field, charge density and potential. Solutions can be obtained in the 
following way: choose as the values of the components v,, v, and v, of the 
velocity v definite functions of x, y and z but containing arbitrary constant 
parameters. ‘This will define a certain type of motion which will be a possible 
one if the arbitrary parameters can be adjusted so as to satisfy the conditions 
given below. Having thus chosen v,, v, and v., form the ‘acceleration 
vector”? A given by 


A=l(v gradv,)+J(v gradv,)+K(v gradv,). (10) 
Then equation (9) will be satisfied provided 
Curl AteeO i” Se ae (11) 


and equation (7) will be satisfied provided 

div (v div A) =0, 
while equations (6) and (8) then serve merely for calculating the field and charge 
density. 

A has to satisfy these two conditions (11) and (12). Ifthe arbitrary parameters 
can be chosen suitably, the resulting values of v and A define a possible space- 
charge motion. 

In terms of v and A other important quantities of the motion are given by 

p=(—me je) divA,-= ~~ ...iv (13) 
E=(-—m/e)A, 
V =(m/e)JAds, 
i=(—me,/e)v div A, 
where iis the current density vector. 


Trajectories, total current and transit times are then also easily calculable, 
as will be seen in examples given below. 


a 
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§3. SOLUTIONS OBTAINED 
Vie Ly) PU IO PIPES Waters he 8 (17) 
DN Ooo Bi) es se (18) 


It is easily verified that equations (11) and (12) are satisfied. Quantities 
of the motion are obtained from equations (13), (14), (15) and (16): 


Example 1. 


Pratt lol Wley/€)imet arses, SIRE Se ty: (19) 
E=—atmje)(lerJy).. 2 1 eines (20) 
V =a*(m/2e)(x? + y*)+constant, © 9 <-> sees (21) 
ba 2a (meen fy). os ee SRL (22) 


It follows from equation (17) that the trajectories are given by the equation 
dy/dx = —4/x,i.e. xy=constant. Thus the trajectories are rectangular hyperbolae 
in planes perpendicular to the z axis and asymptotic to the diametral planes x=0 


and y=0. 


) 


Collector 


R 


———»--— Lines of Electron Flow 
Figure 1. The motion v=a(Ix—Jy). 


A cross section of a system in which such a flow will take place is shown in 
Figure 1. It consists of an infinitely long conducting cylinder of radius R and 
potential Vp, the axis of which coincides with the 2 axis. The cylinder is divided 
into four quadrants by cuts parallel to the axis; two opposite quadrants constitute 
emitters of electrons and the remaining two collectors. The potential on the 
axis Vo, follows from equation (21) : 


oan niet = ee ae, (23) 
This former equation also shows that the equipotentials are concentric cylinders. 


The total current J,,, per unit length of the system is obtained by making use 
of equation (22) to integrate the current flow over the area of unit length of the 


two emitters : Tigp=4@2(megle)R?2. nw ee . (24) 


From equations (23) and (24) follow the relations between current and potential 


difference : TppeSeg(2emy(U RV paVois = wa wees (25) 
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The transit time 7 along a trajectory xy= +47)? (where 7) represents the 
distance of closest approach of the trajectory to the axis) is given by 
T=(Lid)coshe? (Kit). 2 eee (26) 

The interesting result is thus obtained of a flow with transit times varying 
from zero (for 7) = R) to infinity (for 77=0). Finally, it may be noticed, as is seen 
in equation (19), that the charge density is everywhere constant. 

This example has been presented primarily for its theoretical interest. To 
achieve such a system in practice would be difficult, chiefly because of the awkward 
directions of entry of the electrons at the emitters. But an approximation to 
it involving the use of grids as ‘‘emitters”’ with electrons entering from outside 
the system might find applications in space-charge focusing (constant charge 
density), wide-band ultra-high-frequency oscillators (large range of transit times) 
and 90° beam deflections without aid of a magnetic field. 


Example 2. 

Results are presented without detailed derivation; they are deduced in the 
same way as in the last section. Indeed, by a transformation of axes, Example 1 
can be shown to be a special case of Example 2. 


v= l(ay)s-J(O2);., ~~) gay) eaeeeeee (27) 
Asob(leet ly)  «-¢ Sek a ee (28) 
p= —Zablmegieh~ .~ 9 Vs) = 8) Pee (29) 
E= —ab(m/e)(Ix+]Jy), Tae dee (30) 
V =(ab/2)(m/e)(x?+y")+constant, = ...... (31) 

i= —2ab(meoje){I(ay)+]J(bx)}. = 8 — 2. (32) 


The trajectories are given by bx*—ay?=constant and are, therefore, hyperbolae 
since a and 5 must be of same sign, as is shown by (29). 

A cross-section of a system capable of producing such a flow is shown in 
Figure 2. 


Collector 
R 


——»- — Lines of Electron Flow 
Figure 2. The motion v=I(ay)-+J(dx). 
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It is clear that the trajectories are now asymptotic to the diametral planes 


x/b—yr/a=0; BA OE VA) GUL sean (33) 


The potential on the axis V, is given by V)=V,—ab(m/2e)R%, where Vp 
and R represent the potential and radius of the electrodes. The equipotentials 
are again concentric cylinders. 

The total current J,,, is given by J,,,=2ab(a +b)(me,/e)R2. The total current 
is proportional to the three-halves power of potential difference Vp—V,, but the 
relation cannot be put in as perspicuous a form as that of equation (25). 

For a trajectory bx?—ay?=c (c positive) passing through the y=0 diametral 
plane the transit time 7 is given by 


7 =2{4/(ab)}-' cosh“ {4/[b(aR?2+c)/c(a+b)]}}, «see (34) 
and for a trajectory bx? — ay? = —c (c positive), passing through the x=0 diametral 
plane (Figure 2), by 

t =2{4/(ab)}-1 cosh {4/[a(bR?+c)/c(a+b)]}. ws (35) 


Comments on possible applications of this system are similar to those given 
in connection with Example 1. 


§4. DISCUSSION OF METHOD OF SOLUTION 


Consider the case when the velocity components v,, v, and v, may be 
expressed as integral polynomial functions of x, y and z; the dummy suffix 
and Kronecker delta notation are used. Thus, when a suffix appears twice in a 
nomial, that nomial must be summed for all possible values of the suffix, while 
the symbol, 5%” is defined as follows: 84" —1 if l=p, m=q, n=r, and sinn 0 
in all other cases. Further, write S,,,, =x'y™s". Then, for the polynomial case 


Vi La unin t+ JonnnSimn te Kein Simm PYRE nrg (36) 


where the Qn, Opn ANd Cin», ave arbitrary coefficients and /, m and are arbitrary 
indices of the polynomials, which must be so determined that the motion is a 
possible one. [It is to be noted, in accordance with the dummy suffix notation, 
that an expression such as @,,Sjmn, does not mean merely the single nomial 
Ginn Xv", but the sum of terms like it, for all possible values of /, m and n, viz. 
nmr ys”. The values of J, m and n need not necessarily be integral.] 


Hence, from equation (10), 


A = De ln Oia t: Kenn.) 
where 77,,,, iS given by 
Blinc ar Simin Stmanl(E/®) Qin’ a (J) Bym'n’ T (2/2) Crmn’| . 


The two fundamental conditions expressed by equations (11) and (12) can 
now be written as follows: 


Condition 1. 
0/0X(Binn Tian) = 0/0V( Ginn nds 


0/0Y(Cumn Lia) a 0/03(Oimn linen) 
0/02(Gimn Tae) -> 0/0X(Conn 1 s59) 
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Condition 2. 
8/8 Ayan Semin W) + 8/09 Orn Semin W) + 6/02 Crm'n’ Semin: W) = 9, 
where W is given by 
W =0/0%( 41mm Timn) + 9/09 Otmn L mn) +O/02(CrmnL mn): 
Particular results are deduced from the above formulation of the problem. 
in the following manner. 


Case 1. 

Assume er Up Cimn = OS Ding = Oem, 

i.e. assume that the flow is entirely in the direction of the x axis and the velocity 
is a monomial function of the coordinates, since only when J, m and 1 have the 
particular values p, q and r respectively is the coefficient of x'y™z" not zero; 
the application of the two fundamental conditions will show what magnitudes. 
p, q and r must have for a possible flow. 

Condition 1 shows that g=0, r=0, while Condition 2 leads to the equation 
p(2p —1)(3p —2)=Oforp. By referring to equation (13) it is seen that the solutions. 
p=0, p= both give possible motions involving zero charge density and, therefore, 
of no present interest. 

The solution p=? gives the Child—Langmuir flow between infinite parallel 
planes and, for purposes of reference, we give here the characteristics of the flow, 
deducible from equations (13) to (16): 


v=lax'; A=I(2a7/3) x4; 
p=(—2a?/9)(meo/e)x-?; =E=I(—2a?/3)(m/e)x*; 
V =(a?/2)(m/e)x* ; i =1(2a?/9)(me,/e). 
Case 2. 
Assume Bae Ale yyy = adh Dinan = Od pir ys 


i.e. that there is no component of flow in the direction of the z axis, while the other _ 
two components are monomial functions of the coordinates. 

By applying Conditions 1 and 2 it is found that, essentially, only two solutions 
exist involving non-zero charge density, namely, 


p=1, g=0; a 
DD =, ) G =, iY 
a=—b, 
which leads to the motion studied in §3, Example 1, and 
p=09, g= 1, r=(), 
jae eA) r=\(), 
a, b of same sign 


which leads to the motion studied in §3, Example 2. 

It is seen that these two examples represent two of the simplest cases of 
a limited (integral polynomial) class of solutions. By considering binomials and 
other polynomials, and the other classes of function (e.g. fractional polynomials, 
trigonometric functions, hyperbolic functions, etc.), a vast field of knowledge 
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can be opened up. It is to be noted in passing that solutions discarded because 
they involve zero charge-density are still valid solutions and may be of considerable 
interest in low-current or high-voltage phenomena. 

Finally, it must be pointed out that there is considerable scope for pure 
mathematicians in building up existence and uniqueness theorems, establishing 
criteria for solutions with given boundary conditions, studying effects of co- 
ordinate transformations and, in fact, building up a body of theory for electron 
flow comparable to that of potential theory. 
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ABSTRACT. A review of proposals made for altering the formula due to Kénig, on 
which are based measurements of air-particle velocity by Rayleigh disc, leads to the conclusion 
that there is insufficient evidence for altering the numerical constant in the formula. Com- 
parisons between Rayleigh disc and other methods of calibration of individual microphones 
are discussed, and the possibility of difference due to reaction on the sound pressure by the 
vibration of the diaphragm is considered. 


$1. DIRECT EVIDENCE (MEASUREMENTS OF PARTICLE VELOCITY) 

HE familiar formula established by Konig (1891) from hydrodynamical 

theory is now generally accepted as the basis for measurements of sound 

particle velocities, and hence of sound pressures, by the Rayleigh disc. 
The formula is calculated for an ellipsoid in an infinitely extended, uniformly 
moving, non-viscous fluid, and is therefore open to corrections arising from 
practical departures from these ideal conditions. The importance of the 
Rayleigh disc to acoustical measurements has led several investigators to determine, 
by experiment, the magnitude of the small corrections to be applied for such 
of the departures from the ideal as are inevitably encountered in the use of the 
Rayleigh disc for calibrating microphones. In an account of the most recent 
of these investigations Scott (1945) gives a valuable summary cf information 
concerning these corrections—to which his own investigation contributes very 
substantially; but for a correction of the numerical constant in the formula 
proposed by Merrington and Oatley (1939), and to some extent supported by 
Scott, no explanation has been given. ‘The main purpose of the present article 
is to examine the justification for applying this correction. 

K6nig’s formula, as applied to very thin circular discs, states that 


ewe pGeV SUC ee ee ee eee (1) 
where L is the torque on a disc of radius a, set in a stream of fluid of density p 
and velocity V, so that the normal to the disc makes angle @ with the direction 


of flow. 
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Wood (1935) has drawn attention to the fact that in Koénig’s formula the 
velocity V is the relative velocity between the fluid and the disc, whereas with a 
freely suspended disc (as ordinarily used) an alternating movement of the fluid 
imparts a movement to the disc. With sufficient accuracy for the present purpose, 
the effect can be allowed for by introducing the factor (1 — 8a%p/3M) on the right- 
hand side of equation (1), provided that the quantity 8a°p/3M, where M is the 
mass of the disc, is small as compared with unity. This correction has a sound 
theoretical basis (Wood 1935, King 1935), and has received substantial experi- 
mental verification by Wood from experiments made in water, for which, due 
to the greater density, the effect is very much more pronounced than in air. 

Merrington and Oatley (1939) proposed a modification of the formula, which 
not only includes a correction factor for this effect but also alters the numerical 
constant in the formula from 3 to 1:47 (an increase of about 10%). No 
theoretical justification is made for this change, and it is based solely on the results 
of their experimental work. No correction approaching this magnitude was 
found by earlier investigators (Zernov 1908, Barnes and West 1927), who had 
used essentially the same technique, namely, applying the air particle velocity 
by vibrating, at low frequency, a box surrounding the disc and by measuring the 
amplitude of vibration of the box by microscope. 

Zernov has given a reasonably full and complete account ot his experimental 
method and results, as far as his tests with ellipsoids are concerned. He concludes 
from his experiments that Kénig’s formula is fully confirmed for ellipsoids in 
air to within an experimental margin of error of 0:5°%. In arriving at his result, 
however, Zernov introduced a purely empirical correction, amounting to 2:5°%, 
for the effect of the proximity of the walls of the (cylindrical) vibrating box to the 
ellipsoid. He quotes a formula for this correction, but gives no information 
of the data on which he based the formula. The value of correction (2:5%) 
which he used is in the same sense, but is about twice as large as that found by 
Scott (1945, Figure 3). Zernov extended his investigation to discs and found 
discrepancies (due to finite thickness of disc) for which he proposed an empirical 
correction, without quoting the data on which it was founded. This correction 
is also in the same sense, but greater than that found by Scott (1945, Figure 4); 
it is, however, in the opposite sense to that calculated by Konig for a non- 
viscous fluid. The ellipsoids and discs used by Zernov were massive (on modern 
standards, and Wood’s correction mentioned above would be quite negligible. 

The absence of data given by Zernov for his empirical corrections should 
not be allowed to detract from the results of the work on which most of his article 
is concentrated, namely, the experimental verification cf K6nig’s formula for 
ellipsoids. Even if the empirical correction which he made for the proximity 
of the walls of the vibrating box is completely ignored, his experimental results 
lie within 3°% of the value given by theory. From his experiments with discs 
he deduced that, for very thin discs, any correction for alteration of shape from 
ellipsoid to flat disc is very small, that the correction is increased by increasing 
the thickness of the disc, and that it has the opposite sense to that calculated by 
K6nig for non-viscous fluids ; these deductions have received general confirmation 
from subsequent investigation (Scott 1945). 

From their experiments Barnes and West (1927) concluded that, “ provided 
the conditions of use conform within reasonable limits with those on which 
K6nig’s formula is based, then the measurement of acoustic velocities, at any rate 


ee Ah ee ae 
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‘down to the order of 0-1 cm/sec., can be effected by means of a Rayleigh disc 
to an accuracy within 2%” (which corresponds to accuracy within 4% for the 
torque on the disc). Merrington and Oatley discarded the evidence from the 
experiments by Barnes and West on the grounds that two omissions were made 
in their calculations. One-was the correction factor (1—C,) named above, 
where C,= 8a°p/3M. The other was the omission of the increment to inertia 
due to the air in the calculation of the moment of torsion of the suspension fibre ; 
it can also (for convenience of this discussion) be resolved into a correction factor 
(1—C,) on the right-hand side of equation (1) where C, is the ratio J, : J of the 
moment of inertia due to the air to that due to the disc alone. Now J, =16pa*/45 
(Lamb 1906) and J = Ma?/4, whence C, = 64pa?/45M. 

Now information from which the values of C, and C, can be calculated is 
given by Barnes and West for two of the discs used by them for their fundamental 
-calibration, as follows : 


M (gm.) a (cm.) Gi Cy 
First disc 0-:0764 0-6115 0-0096 0-005 
Second disc 0-0622 0-473 0-0055 0-003 


Thus these two corrections, even when combined, amount to appreciably less 
than 2%. 

Even if all evidence other than that presented by Merrington and Oatley is 
discarded it is doubtful whether their experiments were sufficiently complete 
to justify an empirical correction of 10°% to the theoretical formula. There 
is one aspect of their measurement of the air-particle velocity which does not 
seem to have received the attention it requires, and that is the wave-form of the 
‘displacement of the box (and hence of the air) surrounding the disc. The 
measurement of the quantity V?in Konig’s formula was derived from measurements 
of frequency and of maximum amplitude of the displacement of the box. Now the 
presence of a small harmonic component in the wave-form can not only add to the 
mean-square value V?, but it can also alter the maximum amplitude of the wave. 

Take, for example, the case of a sine-wave with maximum amplitude of 
displacement A and angular frequency w; the wave-form of displacement is 
A sinwt and the mean-square velocity is A4?w?. If a small component of third 
harmonic is added, with such phase as to reduce the maximum amplitude, the 
wave-form of displacement having maximum amplitude A is 


A(1+«) sinwt+Ax sin 3ut 


where x is the third harmonic fraction of the fundamental. The mean-square 
velocity A®(1 +.x)?w?+9A?x*w? is greater than that of the pure sine-wave in the 
ratio (1+x)?+9x?:1. Thus, for example, 4% of third harmonic («=0-04) 
can introduce an error in measurement of V? of nearly 10%. Higher crder 
harmonics can, of course, introduce relatively greater errors, for example the 
error with 2° of 11th harmonic could be nearly 9%. 

Now while there may not be sufficient information for asserting that a small 
distortion of wave-form did in fact contribute to the discrepancies between Kénig’s 
formula and Merrington and Oatley’s experiments, there is too much room for 
doubt to permit acceptance of their proposal to alter the numerical constant in 
the formula. In Table II of their paper specific amplitudes of vibration of the 
box are stated, and it is quite consistent that, for tests on any one disc, the dis- 
‘crepancy between experiment and calculation does increase quite markedly 

29-2 
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with increase in amplitude. The mechanism which they used for vibrating the 
box (enclosing the disc) included a valve amplifier and an electromagnetic drive— 
a method which is more open to distortion of wave-form than the simple and 
controlled form of mechanical drive used by Barnes and West. 

At a later date Scott carried out some experiments by a different technique 
which enabled measurements to be made at frequencies up to 4,000 cycles per 
second. He, too, has found discrepancies with Kénig’s formula, though of less 
magnitude than those observed by Merrington and Oatley. His comment that 
the discrepancies which he has observed are too large to be accountable in terms 
of experimental uncertainty is not thought to be sound. It leads to the curious 
conclusion that the numerical constant (3) in Ko6nig’s formula is too small 
by about 7% at low frequencies and by about 4% at high frequencies (2,500 to 
4,000c/s.), but is correct at frequencies between about 1,000 and 2,000c/s.,. 
irrespective of the diameter of the disc. 

Criticism of these two articles is levelled only at the particular conclusion, 
reached in both of them, that the experimental uncertainty is sufficiently small 
to detect and correct for a departure from the theoretical formula, even though 
no theoretical explanation of such a departure can be put forward and established. 
There is much other matter which is of value to the art of acoustical measurements. 
in both articles. 

It is reasonable to conclude that, in the interests of accuracy of measurement,,. 
it is advisable to adopt a technique in which any corrections to be applied to 
Konig’s formula are reduced to very small values. Scott’s paper, in particular, 
contains much useful practical information for guidance on the probable magnitude 
of corrections arising from necessary practical departures from conditions assumed 
in the theory. 

The correction introduced by Wood, 1.e. the factor (1 — C,) mentioned above, 
is a valuable modification of Konig’s formula for application to the calibration 
of a Rayleigh disc. For measurements of sound in air it is convenient to use a 
disc having small moment of inertia, so that the decay factor is not too small and 
the time required for taking readings is not unduly long. A strip of relatively 
massive material, fixed along the vertical diameter of the Rayleigh disc, can 
increase very substantially the mass (and thereby reduce the amount of the 
correction) without appreciably altering either the moment of inertia of the disc 
or the validity of Kénig’s formula. This device is mentioned in §2 of the paper 
by Barnes and West. 


§2. INDIRECT EVIDENCE (MEASUREMENTS OF SOUND, PRESSURE) 


The principal practical application of the Rayleigh disc is to the calibration 
of stable, high-quality microphones for use for the measurement of sound 
pressures. Reference will therefore be made in this section to some comparison 
tests which have been made between calibrations of individual microphones by 
Rayleigh disc and by other methods. Percentage differences which are quoted 
here apply to sensitivities of microphones (that is, the ratio of E.M.F. to sound 
pressure) and not to torque on the Rayleigh disc as in the previous section. 

As a means of practical verification of Kénig’s formula to any high degree 
of accuracy, such tests are less direct and more open to experimental errors than 
those discussed in §1. Uncertainties are introduced, not only by the accuracy 
limitations of the compared method of measuring the sound pressure applied. 
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to the microphone, but also in the accurate measurement of the small £.M.F. 
produced by the microphone, in the dependence on stability of performance by 
the microphone over the interval of time between the two calibrations compared, 
and by possible effects on the sensitivity of the microphone, due to differences 
between the acoustical conditions in which it is used for the calibrations which 
are being compared. (A possible example of an effect of this kind is mentioned 
later.) None the less, the evidence which has been obtained from comparisons 
of calibrations by different methods has relevance and interest. 

Dadson (1946) reports comparisons made by the National Physical Laboratory 
of calibrations of a condenser microphorle by Rayleigh disc and by pistonphone 
oyer the frequency range 60 to 400c/s.; the agreement was within 3%, the 
estimated accuracy of the pistonphone method being to within 3 or 4%. 

Field calibrations of a microphone by pistonphone and by Rayleigh disc 
have been compared in the frequency range 50 to 1,000c/s. by Kennedy and 
Boner (1942). They found substantial agreement at 50c/s. and between 170 
and 1,000c/s., maximum differences being less than 10°% and the average of 
differences practically zero. Between 50 and 170c/s. there were greater differences 
attributed by the authors to ‘‘wobble”’ of the piston. 

For his investigations of effects of obstruction and cavity resonance, Ballantine 
(1930, 1932) measured pressure calibrations by electrostatic actuator, checked 
by thermophone, and field calibrations by Rayleigh disc. At low frequencies, 
where these effects are shown to be small, all the points on his graphs cf the ratio 
of these two measurements lie within 2 or 3°% of unity. 

A direct comparison between the Rayleigh disc and the acoustic radiometer 
methods of measurement of sound energy, made by Hartmann and Mortensen 
(1948), has shown that there is close agreement between these two methods. 

In 1930 a condenser microphone (type D 85880), obtained from the Western 
Electric Company together with the calibration chart of its sensitivity, measured 
by thermophone, was calibrated by Rayleigh disc at the Post Office Engineering 
Research Station. The results of this comparison have not hitherto been 
published. The Rayleigh disc was used to measure both the field sensitivity 
and the pressure sensitivity of the microphone.* The Figure shows the results 
of the two pressure calibrations, by thermophone as supplied by the makers and 
by Rayleigh disc as measured after delivery. It happens that at the lower fre- 
quencies, below about 1,500 c/s., the two curves coincide ; that for the thermophone 
test is level at 0-0027 v/dyne/cm?, and seven points in this range measured by 
Rayleigh disc all lie between 0-00265 and 0-:00275. At higher frequencies the 
Rayleigh disc curve tends to depart by up to about 15°% from the thermophone 
curve, and it appears to reach a maximum sensitivity at a somewhat lower 
frequency. 

It may well be that the differences lie within the combined experimental 
errors of the two tests, plus a possible change in the performance of the microphone 


* The methods and apparatus used for the calibrations are described by West (1932). The 
method used for measurement of field sensitivity has since been changed at this Research Station; 
the pressure sensitivity of a probe microphone (West 1934), which has a reasonably uniform 
response-frequency characteristic, is first measured by Rayleigh disc in the resonant tube; this 
microphone is then used to measure the sound field from a loudspeaker in a relatively large non- 
reverberant room. (As an obstruction to the sound field the probe is no larger, and can well be 
smaller, than the Rayleigh disc and it is much more convenient to use in the open.) ‘The micro- 
phone under test is then measured at the point where the sound field has been calibrated by the 
probe microphone. 
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between the tests. There is, however, also the possibility of a slight difference 
in shape of the frequency characteristics caused by the stiffness reactance imposed 
on the diaphragm by the volume of about 9cm* of hydrogen at atmospheric 
pressure in the thermophone chamber. ‘This stiffness adds to that of the diaphragm 
and to that due to the enclosed air-space behind the diaphragm, and it tends to 
raise the frequency of maximum response associated with the highly-damped 
mechanical resonance of the diaphragm. ‘The argument is speculative in the 
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absence of quantitative knowledge of the mechanical impedance of the diaphragm 
and of the air space behind. ‘This type of condenser microphone has been in use 
in this country and elsewhere as a component of a standard reference system 
for telephone transmission tests; repeated calibration tests carried out over a 
period of many years have shown that individual microphones have a stability 
which is of a very high order. 


83) POSSIBLE EFRECT: ON CALIBRATION, OF DHE ACOUSTICA 
IMPEDANCE TO WHICH THE DIAPHRAGM IS EXPOSED 


As a general comment on the calibration of microphones it is of interest to 
draw attention to the possibility of difference between results for different kinds 
of pressure calibration (particularly at frequencies where the response is influenced 
by mechanical resonance of the diaphragm) due to differences in the acoustical 
impedances to which the diaphragm is exposed in the different measurements, 
if the reaction of the vibration of the diaphragm on the acoustical impedance 
is ignored. 

With a true pressure calibration—where measurement is made of the pressure 
existing at the surface of the diaphragm—the acoustical impedance in front of 


the diaphragm exerts no influence. ‘The author felt the need for experimental. 


% | 


The Accuracy of Measurements by Rayleigh Disc 443. 


verification of this fact, and some calibrations were made at frequencies near 
resonance of a telephone receiver, used alternatively as a microphone and as a 
source of sound, with different acoustical impedances imposed by forming an 
enclosure in front of the telephone receiver and by altering the volume of this 
enclosure. A calibrated probe microphone was used for measuring the sound 
pressures at a point just in front of the diaphragm. A source of sound formed 
part of the wall of the enclosure for generating sound pressures therein when the 
telephone receiver was being calibrated as a microphone. Electrical measure- 
ments were made in terms of voltages at the receiver terminals. The mechanical 
resonance of the diaphragm showed, as expected, a very pronounced peak in all 
the tests, but whereas the frequency of maximum sensitivity changed from about 
1,300 to about 1,000 c/s. for variations of enclosure volume from 1-5 to 20cm? 
when the telephone receiver was calibrated as a source of sound, it remained 
constant at about 1,000c/s. when the telephone receiver was calibrated as a 
microphone. (Precautions were taken to eliminate changes due to temperature 
variations.) 

For the microphone calibration, the pressure iz the enclosure may be regarded 
as made up of two components, one causing the diaphragm to vibrate and the other 
resulting from its vibration. The relative magnitudes and phases of these 
components alter as the frequency passes through that of resonance of the 
diaphragm. It would appear that the analogue, using the impedance concept 
applied to electrical circuits, is that of a generator, having an internal reactive 
impedance, applied to a series-tuned circuit (the diaphragm); the response of 
the tuned circuit, in terms of the ratio of current flowing in it to the potential 
difference across its terminals (i.e. ratio of velocity of diaphragm to’ pressure in 
enclosure), is independent of the impedance of the generator. ‘The two compon- 
ents of sound pressure in the enclosure are represented, respectively, by the 
E.M.F. and the potential difference across the internal impedance of the generator. 
The corresponding analogue in the case of the calibration as a source of sound 
would be to regard the E.M.F. as generated in the tuned circuit (which has the 
reactive impedance of the enclosure in series with it) and the response of the 
mechanical system as the ratio of the potential difference across the impedance 
of the enclosure to the E.M.F. (i.e. of pressure in enclosure to force on diaphragm). 

Some of the methods used tor measuring pressure sensitivities of microphones 
do not yield true pressure calibrations (in the sense used here). Thus, for example, 
in the thermophone and the electrostatic actuator methods it is usual to ignore 
the contribution to the pressure on the diaphragm attributable to the movement 
of the diaphragm, and, therefore, to include the effect of the acoustical impedance 
to which the diaphragm is exposed in the calibration. It is true that at fre- 
quencies remote from that of resonance of the diaphragm (or in cases where the 
resonance is so highly damped as to exert no appreciable influence on the response) 
the contribution to the pressure due to the movement of the diaphragm is very 
small; in other words the load imposed by the acoustical impedance is an insigni- 
ficant fraction of the mechanical impedance of the diaphragm. In fact the 
effect has been considered and quantitative estimates of its magnitude made by 
Ballantine (1932) for the thermophone method of calibration, but consideration 
was limited to frequencies at which the diaphragm is stiffness-controlled: it did 
not include the increase of magnitude and change of phase which occurs at 
resonance of the diaphragm. 
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The Rayleigh disc method of calibration in a resonant tube does yield a true 
pressure calibration, since any contribution to the stationary wave in the tube 
caused by vibration of the diaphragm is included in the measurement made by 
the Rayleigh disc. Thus, if mechanical resonance of the diaphragm does exert 
appreciable influence on the response of a microphone in the range of frequencies 
tested, it is to be expected that the maximum response due to the resonance will 
occur at a somewhat higher frequency for the thermophone than for the Rayleigh 
disc calibration. Moreover, the sensitivity at the point of maximum response 
depends largely on the amount cf mechanical resistance opposing the movement 
of the diaphragm (in the case of a condenser microphone the resistance resides 
mostly in the small air space between the electrodes). Since the resistance 
tends to increase as the frequency is raised, a somewhat smaller maximum response 
is to be expected when the maximum occurs at a higher frequency. The result 
is a difference (though possibly only slight) in the shape of the frequency-response 
characteristic as measured by these two methods of calibration. 

In the case of a measurement of field sensitivity, the sound pressure measured 
is that existing in the free air at the position of, but in the absence of, the micro- 
phone. There is thus no question of movement of the diaphragm contributing 
to the value of the pressure measured, and the acoustical impedance load on the 
diaphragm is that due to the free air—mostly a small inertia reactance. 
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ABSTRACT. A new type of gravimeter is described, which depends for its operation on 
the change of frequency of the natural vibration of a vertical wire stretched by a weight. 
The advantages of such an instrument are presented, and a description given of the design 
which has so far been evolved. The results of laboratory tests and of trials carried out at 
‘sea in H.M. Submarine Talent are given. 


§1. INTRODUCTION 


T present the pendulum is the “standard” instrument for relative 
measurements of gravity: all other gravimeters have to be checked against 
this standard more or less frequently. This is usually effected by 

visiting known gravity stations at intervals, and working thereafter in a closed 
loop each day, so as to allow for the daily drift of the instrument. 

One of the aims in making the vibration type of gravimeter was to provide 
a substandard to the pendulums. The cause of drift in nearly all other existing 
gravimeters lies in the spring: the extension of the spring is directly proportional 
to the total attraction of gravity. Since the elastic properties of most materials 
vary slowly with time, the zero reading of the gravimeter will alter, and this limits 
the time over which the instrument may be used to obtain results accurate to 
any given amount. It is possible, however, to design a gravimeter in which the 
elastic properties of the materials do not affect the readings directly. Such an 
instrument has been described by Bertrand (1938). It consists of a mass 
suspended from a fixed support by a thin wire. If the wire is made to vibrate 
transversely, its natural period will depend on the tension, and hence on the value 
of gravity. 

The frequency of the real wire differs from that of an ideal, perfectly flexible, 
inextensible wire for three reasons: firstly, the wire has a finite rigidity; secondly, 
the length of the wire varies with the tension; and thirdly, the supports at each 
end of the wire are not absolutely rigid. Suitable design renders these terms 
small, and drift due to variations in them can be largely eliminated. ‘The fre- 
‘quency of the ideal string depends inversely on its length and the square root 
of its mass per unit length and directly on the square root of the tension. The 
tension in the string depends on the mass suspended by it, and this and the mass 
of the string can easily be kept constant to 1 in 10’ provided precautions are taken 
to avoid corrosion. The length of the string, apart from the effect of creep, 
depends only on the temperature, and the apparatus is accordingly placed in a 
thermostat. Experiment shows that the effect of creep can be reduced to less 
than 1 part in 10° per day. The frequency is also dependent on the air pressure, 
as in the case of pendulums, and the apparatus is therefore evacuated. 

The principal advantage which the instrument possesses over an ordinary 
gravimeter is its extended scale. ‘The time needed for an observation to attain 

* Patents pending. 
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a desired degree of accuracy depends only on the accuracy with which the phase 
of the wire can be measured with respect to a fixed frequency standard. At 
present this is about 1/10 cycle and, since the string vibrates at 1,000 c/s., a reading 
correct to 1 in 108 can be obtained in 100seconds. This accuracy is unaffected 
by the actual frequency difference between the wire and the standard; thus 
variations in g of 2%, such as are encountered over the earth, can be measured 
with the same accuracy. The instrument is therefore suitable for geodetic use. 

The accuracy of phase-comparison cculd easily be improved to 1/50 cycle 
and thus an accuracy of reading of 1 in 108 obtained in just over half an hour.. 
The instrument could then be used to measure earth tides, provided that its 
drift were uniform. 

The instrument also possesses integrating properties. Since the beat- 
frequency is, to the first order, proportional to g, then the total number of beats 
measured in any given interval is proportional to the mean value of g, whatever 
may have been the form cf its variation during the interval. This property makes 
the vibration gravimeter especially valuable for measurements in submarines. 
Indeed the only device which shares this property is a pendulum-apparatus, 
which so far has formed the only successful method of measuring gravity at sea. 

The instrument is unlikely to compete with static gravimeters for ordinary 
commercial surveys, since its associated electronics and power supplies render 
it inconvenient. Nevertheless, it might be possible to make a model small 
enough to be lowered down a borehole. A well-known strain gauge depends 


on the principle outlined, but is not nearly so accurate as the instrument here 
described. 


§2. DESIGN 


In his original paper, Bertrand described the apparatus as follows: 

‘‘A vertical string is tensioned by the weight Mg of a mass M. The fre- 
quency of its fundamental F=(1/2L),/(Mg/m) (L is the length of the string, m 
its mass per unit length) varies with g and enables variations therein to be measured. _ 
From the measured frequency one can deduce the value of g, using a constant 
applicable to the apparatus, which constant can be found by adding known weights 
to M.” 

A schematic diagram of the instrument is shown in Figure 1. The weight 
is suspended from a beryllium copper strip about 5cm. long. This particular 
instrument was designed for use at sea, and it was therefore desirable that the 
weight hanging on the string should be damped, so as to prevent any pendulous 
swinging, with consequent inaccuracy in the determined value of gravity. The 
weight was made in the form of a copper disc, with a downward projecting rim, 
and a boss in the centre in which the half-cylinder at the end of the string could 
be clamped. A damping magnet was made, with four sets of poles symmetrically 
placed round the weight, which produced four horizontal radial fields, cutting 
the rim of the weight, and four vertical fields, cutting the disc portion of the 
weight. ‘The latter was thus damped for lateral movements, twists about the 
axis of the disc, twists about any radius of the disc, and vertical movements. 
The constancy of the mass of the weight is assured, for copper does not readily 
corrode and has a low vapour pressure. 

The string must have a cross-sectional area sufficient to support the weight ; 
a tension of approximately one-fifth of the breaking stress is used. To the first 


A Dynamic Gravimeter of Novel Design 447° 


order, the frequency of a flexible bar, of finite rigidity, clamped at each end, is 
given by Rayleigh (1929) as 
n=7al/I[1 +2(b/a)(«/L)], 
| so long as 2(b/a)(«/I) is small. : 
Here a= /(T/p); Tw is the tension in the string, pw the linear density, 
w the area of cross section of the string, / the length, b= /[(q¢+T)/p]; ¢ is the 
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Figure 1. Schematic diagram of gravimeter (clamps omitted). 


Young’s modulus of the string, and « is the radius of gyration of the section of the 
string about its median plane. n=va/l is, of course, the frequency of the ideal 
string. 

In order that the frequency shall depend but little upon the elastic constants 
of the wire 2b«/al must be small compared with 1. Now b/a=,/[(q+T)/T], 
and T cannot be increased indefinitely, or the wire will creep, and eventually 
will break. A safe limit is about T/q=1/2,000. 
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« may be made small by using a thin flat wire since x =d/24/3 where d is the 
thickness of the wire. 

In the practical case, we have g= 1-27 x 10! dynes/cm?, Tw =65 x 981 dynes, 
.w=125 x 10-§cem?, d=5 x 10-3cm.; then «= 1-44 x 10-°cm., /=5-2cm. Hence 
T =65 x 981/125 x 106 dynes/cm?, and \/[(¢ + T)/T] =50, giving 2bx/al = 2-8 x 10-?. 
Thus the rigidity correction is only 2-8°, of the total frequency, and the formula 
may legitimately be used. 

It is clear from the above that the wire should be flat, and owing to the ease 
with which it may be hard-soldered, beryllium copper was used. ‘The available 
‘suitable strip was 0-010in. x 0-002in. in cross section: using this size with a 
string approximately 5 cm. long, and a weight of approximately 65 gm., a resonant 
frequency of nearly 1,000c/s. was obtained. The wire hangs vertically between 
the poles of a permanent magnet, the lines of force being horizontal and parallel 
to the plane of the strip. If the wire vibrates across the lines of force, an alternating 
E.M.F. is generated between its ends: conversely, if an alternating E.M.F. is 
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Figure 2. Circuit diagram of driving circuit for gravimeter. 


‘applied, the amplitude and phase of the vibrations, and hence the current in the 
wire, depend on the frequency. ‘The wire accordingly behaves as a resonant 
circuit element, but of rather low dynamic impedance. ‘The wire is attached 
at the top to a block, which is adjustable for tilt and is insulated from the main 
body of the instrument; one connection is made to this block, and the weight is 
connected to another terminal by means of a fine beryllium copper hair-spring. 

The circuit of the oscillator used for driving the string is shown in Figure 2. 
‘The voltage appearing across the string is amplified, and the amplified voltage 
fed back to the string. In order that the circuit shall act as an oscillator, the 
energy fed into the wire must be equal to the energy lost : owing to the low dynamic 
impedance of the string, this renders a two-valve amplifier necessary. The ampli- 
fier will ‘‘pull” the frequency of the string slightly from its natural value; the 
‘Q of the string is of the order of 20,000, and the phase change through the amplifier 
as about 1°. Thus the total pulling of the frequency is only 1 in 108, and variations 
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in this will be less than 1 in 107. The pulling of the frequency makes it undesirable 
to use a transformer in the circuit, for if any direct current in the transformer 
should vary, the inductance will alter, and with it the phase change. Thus an 
extraneous frequency change would be observed. 

The output frequency was compared with a crystal frequency standard.. 
The method of comparison was suggested by members of the G.P.O. Research 
Laboratory, Dollis Hill. The thousand cycle signal from the crystal clock is. 
split into three phases by condensers and resistances. Each phase is made to 
produce beats with the thousand cycles from the oscillator, thus producing three 
phase beats, which are applied to a small ‘“Selsyn”’ motor. ‘The motor makes 
one turn for each beat, and fractions of a beat may easily be estimated. ‘The 
beats are counted by a gear train and dials. This method has the great advantage 
over most others, that it differentiates between frequencies greater and frequencies 
less than the comparison frequency: in the one case the Selsyn rotates clockwise, 
and in the other, anticlockwise. This is a very valuable feature for sea use, where 
vertical and horizontal accelerations may vary the apparent value of gravity 
by nearly one per cent. For, if the apparatus did not discriminate between fre- 
quencies greater or less than 1,000c/s., then the nominal beat frequency would 
have to be greater than the maximum variations caused by the movement of the 
ship. 

The integrating property of the instrument was pointed out above. If the 
apparent value of gravity in the submarine is fluctuating then the mean value 
measured, g, is to the first order 


dx/dt) dx/d 
7. gdt=g+ a t)o 


The first-order correction arising out of the difference of vertical velocity 
(dx/dt), —(dx/dt), at the beginning and end of an observation could be eliminated 
if the dials were always photographed when (dx/dt),—(dx/dt),=0. This could 
be done by obtaining first a voltage proportional to d*x/dt® and integrating: 
the camera could be arranged to trigger automatically whenever (dx/dt), = (dx/dt),. 
The technical difficulties of doing this were not solved in the present instance, 
however. The second-order correction remains in any case and is discussed 
below. 

An electronic thermostat was constructed. An a.c. voltage was applied 
to a Wheatstone bridge, the arms of which consisted of coils, of different tempera- 
ture coefficient but equal resistance at the balance point. The output from the 
bridge was amplified and applied to a phase sensitive rectifier, which discriminated 
between too high and too low a temperature, and controlled the heater accordingly.. 
The coils were wound on a large copper can, two of 47s.w.G. copper wire, and 
two of 32s.w.G. constantan wire, thus providing four coils of the same number 
of turns, and with approximately the same resistance, but with temperature 
coefficients 0:0043°c. and zero respectively. The coils were wound similarly 
so as to balance out their self-inductance, and the bridge was wound so that there 
was no mutual inductance between the input and the output points. ‘The heater 
was also wound non-inductively, to avoid stray magnetic fields. The 1,000c/s. 
voltage from the crystal frequency standard was used for the a.c. supply, since 
it had a good waveform; a signal of 5 volts was applied to the bridge. ‘The 
out-of-balance voltage for an increase in temperature of 0-001° c.. from balance 
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was 21lpy.; this was amplified by a factor of about 5 x 104, and when rectified, 
gave about 1 volt across a polarized relay, which was more than adequate to 
operate it. In practice, the hunting of the temperature was estimated at about 
-0-0005° c. on either side of the balance point. 

The strip of beryllium copper was originally clamped at each end between — 
blocks of brass, as in Figure 3, but no very accurate provision was made for 
ensuring that the edges of the blocks clamping the wire were 
level. This meant that, for half the cycle, the wire was 
effectively of a different length from that of the other half- 
cycle, which may be more clearly seen in the exaggerated 
sketch. Alternative methods of fixing were investigated. 
The next method used was to soft-solder the wire to two 
small half-cylinders of brass, with the wire lying along the axis 
of the cylinders. These half-cylinders then located the ends 
of the wire, and were clamped between the blocks, a ‘“‘vee”’ 
groove being cut in one of the blocks to locate the half- 
cylinder. As the wire was soldered to the half-cylinders 
whilst held in a jig, this method ensured that the wire hung 
vertically, within the accuracy of machining. Experimental 
results encouraged us to proceed further to eliminate the zero 
drift, which was most pronounced when the weight was 
clamped and unclamped. This led to the possibility that 
when tension was applied suddenly on unclamping, the wire 
might slip slightly through the soft solder. Accordingly 
attempts were made to hard-solder the wire; this was 
finally accomplished in an atmosphere of nitrogen in an electric furnace. After 
soldering (in a jig) the wire was quenched, and then heat-treated to re-harden 
the beryllium copper strip, which had been softened by the high temperatures 
involved. 

The zero-drift with this string was indeed much less than that with the soft- 
soldered strings, but still was not negligible. This was felt to be due to the method _ 
of clamping, whereby the weight was lifted, and the tension removed from the 
string. It was therefore arranged to clamp the weight, yet keep the string under 
tension. In order to do this, some part of the assembly was mounted on a spring 
to prevent very slight movements of the weight when clamped from stretching 
the wire. This was accomplished by mounting the top end of the string on a 
spring. Ifsuitably orientated, the spring will present a high mechanical impedance 
to the top cf the wire, and yet allow vertical movement, which is the desired 
feature. Great difficulty was at first experienced in soldering the wire to the top 
end of the spring, owing to oxidation of the beryllium copper, even when in an 
atmosphere cf nitrogen. A hydrogen furnace was used to solder the wire, but 
was later discarded in favour of a very small oxy-hydrogen flame. This enabled 
soldering to be effected without oxidizing the rest of the wire by overheating. 
‘The wire was subsequently heat-treated in the hydrogen furnace. 

A clamping device was considered which could be operated without disturbing 
the apparatus. This consisted of a tube carrying two arms with rounded ends 
which press upwards on the weight when clamped. The tube slides in a vertical 
guiding tube, and as it moves upwards it cperates two levers through a system 
of push-rods, each carrying a ball-bearing which moves downwards toward the 
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weight. Consequently, the weight is clamped in the position in which it normally 

rests ; indeed the frequency of the wire is very nearly unchanged when the weight 
is clamped. ‘The motion of the clamps was produced by a “‘ Sylphon”’ bellows, 
the top surface of which formed part of the wall of the vacuum system. When 
‘the under surface of the bellows was at atmospheric pressure, an upward force 
was produced on the bellows, and transmitted to the clamps. On evacuating 
the under compartment, the clamps were returned to their free position by 
“springs. 

Before the spring at the top of the wire was used, it was noticed that the 
un-sprung wire acted as a vertical seismograph. This was discovered by acci- 
dentally dropping a large weight on the floor. The cause was easily ascertained : 
for the tension in the string would change momentarily, although the mean tension 
over a time would remain constant. ‘The frequency depends on the square root 
of the tension; therefore if the tension varies, the mean frequency will be lower 
than it would be under static conditions, by a correction of the second order of 
the variation of frequency. ‘This is also the explanation of the second order 
correction due to vertical accelerations of a submarine at sea. The spring, 
however, will ensure that the tension in the wire remains more nearly constant 
during sudden vertical accelerations. 


§3. LABORATORY TESTS 


The effect of temperature on the frequency of the wire was found to be 
0-:0049 c/s. per °c., in 1,000c/s. This corresponds to approximately 10 milligals 
per ° C., over a range from 15°.c. to 27° c., and is only about one third the coefficient 
expected from the known coefficient of expansion of the wire and the variation 
of its elastic properties with temperature. 

The pressure effect appears to be due to three causes: firstly, the air has a 
buoyancy effect; secondly, a certain amount of air is carried with the wire, which 
increases the effective mass per unit length; these predominate at the higher 
pressures between 0-1mm. and 760mm.; thirdly, by virtue of its damping pro- 
‘perties, the air lowers the QO of the wire, and hence the frequency changes, due 
to the fact that there is a phase-change in the amplifier. This third cause 
becomes apparent only at low pressures (below 0-1mm.) where the first two 
__ effects are small. The variation of frequency with pressure is shown in Figure 4. 
_ The dynamic impedance of the string increased steadily as the pressure was 
reduced, a graph of the impedance plotted against the logarithm of pressure 
being a straight line within the limits measured. The impedance rises from 
24 ohms at atmospheric pressure, to 80 chms at 2 x 10-?mm. 

The drift experienced when the instrument was in use may be seen from the 
following results, obtained when the instrument was mounted in H.M. Submarine 
Talent. 

Fort Blockhouse 19th May Beat frequency 2-8607c/s. 
af és 21st May :f y 2SO2n C/s: 
Portland 27th May » i 2°866 c/s. 

A change in beat frequency cf 0-001c/s. is equivalent to —2 milligals: the 
value of gravity at Blockhouse is 2 milligals greater than that at Portland. 

The amplitude effect appears to depend on the cube of the amplitude: it is 
easier to maintain the amplitude of oscillation of the wire constant than to apply 
acorrection. This is effected by means of an automatic volume control system 
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in the oscillator, which maintains the voltage across the string at a constant value.. 
In practice, this voltage is about 5 microvolts, which corresponds to an amplitude 
of about 10-5cm. The noise level is of the order of a tenth of a microvolt: but 
the amplitude of oscillation was not reduced further, as the automatic volume 
control was sufficiently good to keep the error due to this cause below half a 
milligal. (The frequency is actually 0-005c/s. lower than the frequency at 
zero amplitude when the apparatus is operating normally.) When the mains 
voltage was varied from 160 to 195, the frequency increased by less than 0-0005 c/s. 
in 1,000c/s. Tilting the apparatus appeared to vary the frequency approximately 
as 1/cos 0 if @ was small, a tilt of 1° from the centre frequency minimum causing 
an increase of frequency of 0-15c/s.in 1,000c/s. So long as the levelling is correct 
to within 0-05°, error due to this is less than 1 milligal. 
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So far, the only field tests made on this instrument have been at sea during 
the recent cruise of H.M. Submarine Talent in the English Channel. A Vening— 
Meinesz 3-pendulum apparatus was used as a standard of reference. The 
instrument was mounted in gimbals to eliminate the effect of the slight rolling 
of the vessel even when submerged. Records were made automatically at 5-second 
intervals with a 16mm. camera, which photographed together the beat-counter 
and a synchronous clock governed by the crystal frequency standard. 

Several difficulties were encountered: electrical interference made it difficult 
to keep the string in oscillation: the apparatus accidentally became filled with 
oil while the pump was stopped, and the wire was damaged while adjustments. 
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were being carried out and had to be replaced. Despite these accidents, several 
readings were obtained which proved to be in very reasonable agreement with the 
pendulum apparatus. 

The errors involved whilst at sea were mainly due to two causes: (a) imper- 
fections in the gimbal system, which caused the gravimeter to be tilted by the 
rolling of the submarine, and (b) the movement imparted to the submarine by 
waves, which may vary the apparent value of gravity by one per cent with a period 
of about six seconds. The maximum estimated error due to the first cause is 
3 milligals, and that due to the second, 1 milligal. The probable error of a 
pendulum observation is 1-5 milligal. The results obtained are shown in the 
table below. 


Station Value of Reading Error Station Value of Reading Error 
number gravity of [mgals.] | number gravity of [mgals. } 
[pendulums] gravimeter [pendulums] gravimeter 
84 981-080 981-076 —4 Cherbourg 981-043 981-042 —1 
85 981-064 981:066 +2 93 981-060 981-055 —5 
86 981-073; 981-072 —1:5 [94] 981-086; 981-080 [—6:5] 
87 981-043 981-040 —3 95 981-085 981:090 + 
88 981-034 981-031 —3 96 981-076 981-078 aP2 
89 981-088 981-090 +2 97 OS 1A Oe SitatOG —4 
90 981-112 981-112 0 52h 981-105 981-108 +3 
91 981-106 981:102 —4 522 981-104 981:108 +4 
92 981-078 981-080 +2 [527] 981-098 981:106 [+8] 


Root mean square value of errors=4. 


The two results given in brackets were taken under unusually rough conditions ; 
if these are discarded the r.m.s. error is 3-2 milligals. The standard error of 
the pendulums is 2-4 milligals, giving a standard error of the same order for the 
gravimeter, and a probable error of 1-5 milligals for the gravimeter alone. 


§5. CONCLUSION 

Although only in the initial stages of development, the apparatus already 
shows great promise. ‘There is still a small drift of frequency with time, but 
there are indications that this grows smaller as the wire ages. The results obtained. 
at sea during the cruise in H.M. Submarine Talent show agreement between the 
pendulums and the gravimeter within the limits of experimental error of both. 
This is a great step forward as hitherto the use of pendulums has been the only 
successful method of measuring gravity at sea. ‘The performance could be 
considerably improved by using better gimbals, and readings simplified by the 
use of further electronic aids. 

If the drift can be reduced sufficiently, the apparatus will be of geodetic use 
(i.e. large variations of gravity may be accurately measured), and may even prove 
more accurate than pendulums. Investigations on the densities or rocks, by 
lowering the instrument down a borehole, should prove practicable. In order 
that this could be done, the apparatus would have to be permanently evacuated, 
the clamping mechanism changed so as to be operated from the surrace, and the 
gimbals improved. The problem of thermal control of the instrument might 
be overcome by good thermal insulation and a recording thermometer. Owing 
to the great precision with which the frequency may be measured, the apparatus 
may have uses for investigation of earth tides at suitable permanent stations. 
The possibility of using it as a seismograph has to be looked into further. 
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ABSTRACT. In order to determine whether any structural change takes place when a 
liquid is supercooled, measurements have been made on the viscosity of pure phenyl ether, 
which by its melting point (26-85° c.) is suitable for accurate work. It has been found that 
in both the ordinary and the supercooled state the liquid obeys a formula of the type 
=A exp (b/T), but that the value of } in the supercooled state, which was investigated down 
to 17° c. below the melting point, is significantly larger than the value of 6 for liquid above 
the melting point. It is shown that an accurate value for the melting point can be obtained 
from viscosity measurements without the liquid ever becoming solid. The temperature 
coefficient of density shows no detectable change in passing through the melting point. 


$1 INTRODUCTION 


HE phenomenon of supercooling is well known. A substance can often 
exist in the liquid state even though its temperature may be several degrees 
below the melting point of the substance, and even though the liquid be 
stirred or made to flow. ‘The liquid may remain supercooled indefinitely provided 
no solid crystal of the substance is dropped into it. The melting point of the 
substance may, however, be quite definite, as shown for example by the cooling 
curve for phenyl ether (Figure 1) taken with frequent additions of tiny crystals of 
the substance. The horizontal part of the curve gives the melting point at 26-85° c. 
It is of interest to know if any change in the nature of the liquid occurs as it is 
cooled and passes through the melting point into the supercooled region. A very 
suitable property for investigation of any such change in the liquid is that of 
viscosity, for on theoretical grounds a pure, non-associated liquid obeys the law 
given by Andrade (1934), viz. y= Aexp(E/RT), where 7 is the viscosity of the 
liquid at absolute temperature 7, Rk is Boltzmann’s constant, and A and E are 
constants for the particular liquid, whilst experimentally, relative viscosities may 
be measured to a very high degree of accuracy of the order of one part in 100,000. 
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S23 AR PARA US 


A pipette viscometer was used, in which a fixed volume of the liquid is drawn 
up into the pipette and allowed to discharge into a reservoir (Figure 2). Any 
change in working volume by evaporation is avoided by the incorporation in the 
‘suction tube of a sphere S containing some of the liquid at the same temperature as 


Temperature (°C) 


0 5 10 15 20 
Time (minutes) 


Figure 1. 


that inthe viscometer. Local distillation is prevented by immersing all glass parts 
having contact with the liquid or its vapour in a thermostatically controlled water 
bath, level LL. All ground joints are made water-resistant by the application of 
silicone grease. With a large, well-stirred water 
bath, the temperature can be maintained constant 
to within 0-005° c. which, for a time of flow of about 
300 sec., is equivalent to a variation in flow time 
of 0-03 sec. Since the extreme deviation to be 
expected, using a stop watch to measure the time 
of flow ,is 0-3 sec., this method was discarded in 
favour of a photoelectric method. The liquid is 
drawn up above the lower bulb and by adjustment 
of the tap is allowed to fall under its own hydro- 
static head. An image of an illuminated slit is 
formed at A and the light passes on through the 
liquid and is focused on to a photoelectric cell 
controlling a high speed relay through a suitable 
amplifier. As the liquid meniscus passes A it 
scatters the light momentarily and the cell operates 
the relay, which starts an electric clock. A similar 
interruption of a beam of light when the meniscus 
reaches B is used to stop the clock, which is of the 
synchronous motor type controlled by a tuning 
fork, and has four dials reading down to 0-01 sec. 
The mean deviation for ten successive times of flow 
never exceeded 0-03 sec., showing that full use was 
being made of the constancy of the temperature 


Figure 2. 


30-2 
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provided by the thermostat. Over a limited range of temperature the working” 
volume changes only slightly, due to thermal expansion; therefore, rather than 
attempt to keep the working volume constant, the change in time of flow at 
a given temperature, produced by adding small volumes from a 1 cm? pipette, 
was found, and an appropriate correction deduced from the known coefficient 
of expansion of the liquid. The correction was never more than 0-5%. At the 
higher temperatures, when the rate of flow is some three times the lowest value, 
a small correction (again never more than 0-5°%) to the observed time of flow was 
applied to allow for the kinetic energy of the issuing liquid. The ratio of the times 
of flow at two temperatures, for the lower of which the correction was certainly 
negligible, was determined for various values of the hydrostatic head. This ratio 
tended towards a constant value as the head diminished and the times of flow 
correspondingly increased. From this limiting value the correction was deter-- 
mined. 
§3. DETERMINATION OF THE DENSITY 


For a liquid flowing under its own hydrostatic head, the viscosity is proportional 
to the product of the time of flow tand the density p. The density was determined 
at various temperatures, for both the normal and the supercooled regions, using a. 
dilatometer. The volume occupied at any temperature by a fixed mass of liquid. 
was determined by cathetometer observations on the position of the meniscus in the: ~ 


6 V obs Veale V obs— V cate 6 V obs Veale Vobs— V cate 
(65) (cm?) (cm?) V (2xek), (cm?) (cm?*) V 
12:40 10:6667 10:6664 -+2-8x10- 25-98 10-7846 10-7849 —2-8x10-> 
13-75 10-6782 10-6780 +1:9x10-5 26°85 Melting Point 
14:07 10-6811 10-6808 +2:8x10-5 28-01 10-8025 10-8029 —3-7x10-> 
16:04 10-6978 10-6978 0-0x10- 30:17 10-8224 10-8222 +1:9x10-5 
13-139 510-7157 10-7159 = 1-941 0 32:24 10-8405 10-8407 —1-9x10-° 
19-82 10-7308 10-7307 +0:9x10-5 33:82 10:8548 10-8551 —2-8x10-° 
21:17 +10°:7427 10-7425 +1-9x10-5 35:27 10:8670 10:8679 —8-3x10-> 
21:98 10°7499 10-7496 +2-8x10-5 36:02 10°8740 10°8748 —7:-4x10-5 
24:12 10°7684 10-7685 —0-9x10-5 


capillary tube, the usual allowance being made for the expansion of the vessel and 
for non-uniformity of the bore of the calibrated capillary tube. In the Table the 
observed volumes of 11-5349 gm. of phenyl ether at temperatures @° c. on each side 


of the melting point are given, together with the error obtained by fitting these 
values to the equation 


V =10-5612(1 +7-9188 x 10-46 +. 9-033 x 10-7 62). 


Over the whole range these errors are small and unsystematic, showing that there is 
no abnormal density change as the liquid passes through its melting point. 


§4. VISCOSITIES 


Since 7 = K(tp), where K is an apparatus constant, then if 7 =A exp (E/RT), a 
graph of log 7 against 1/T should give a straight line. The results for phenyl ether 
(M.Pt. 26°85° c.) are plotted in Figure 3. The points fall into two groups corre- 
sponding to observations made above and below the melting point, the points of 
each group lying on a straight line, but the line for the supercooled liquid being 
steeper than that for the liquid above the melting point as shown in Figure 4. 
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‘Figure 4. Graphs of log y against 1/T for phenyl ether above the melting point (outer scales) and 
supercooled (inner scales). 
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The increase in slope is about 13%. To show that the change of slope is abrupt, 
the equation to the straight tine for the liquid above the melting point was found by 
the method of least squares and the deviations of all experimental values from the 
values given by this equation plotted against 1/7 as shown in the inset to Figure 3. 
Above the melting point the deviations are small and both positive and negative 
about the zero line, but in the supercooled region the deviations are all positive and _ 
increase progressively as the temperature is lowered. The point of intersection 
of the two straight lines shown in Figure 3 is approximately at the melting point, 
but a more rigorous test was applied as follows. ‘The equation to the best line 
through the points corresponding to temperatures above the melting point was. 
calculated by the method of least squares, and a corresponding equation for the 
best line through the points in the supercooled region. By solving these equations 
algebraically the temperature corresponding to the point of intersection of these 
two lines is 26-95° c., which is in good agreement with the value 26:85° c. obtained 
for the melting point by the cooling-curve method. The striking fact emerges that 
the melting point of a substance can be accurately determined from viscosity 
measurements made into the supercooled region without the liquid ever becoming 
solid! In connection with Figure 3 it should be mentioned that in order to avoid 
any systematic error arising in the apparatus, the times of flow were measured. 
first in the region above the melting point, then in the slightly supercooled region,, 
back into the normal region, further into the supercooled region, and so on. 


§5. DISCUSSION 


The fact that the graph of log 7 against 1/T gives one straight line for the heated. 
liquid and another for the supercooled liquid is interesting since it shows that both 
approximately obey the Andrade Jaw and therefore both are normal liquids. ‘The 
change in slope indicates that the value for E in the Andrade equation increases as. 
the liquid passes through the melting point into the supercooled region. ‘The 
magnitude (13°) of the increase in the value of £, which on Andrade’s theory 
represents an energy of interaction of the molecules, is quite striking. For com- 
parison it may be noted that the isomeric change from CH,—O—C3H, to 
C,H;—O—C,H, causes a change in the value of EF, as deduced from viscosity 
measurements, of only 1:3% (see Andrade 1934). An appreciable structural. 
change therefore sets in at the melting point, and it is hoped that further evidence: 
in support of this may be obtained from a comparison of the x-ray diffraction 
patterns produced by the liquid above and below the melting point. A careful. 
examination of the (log, 1/7) curve for the supercooled liquid indicates that the 
curve departs slightly from a linear law at the lowest temperatures, being slightly: 
concave upwards. ‘This curvature seems to indicate a progressive change in 
structure as the liquid is further supercooled. 

‘That the effect is not due to the formation of the well-known liquid crystals. 
(e.g. see Lehmann 1889) was shown by observing the cooling liquid with polarized 
light under a microscope. No restoration of the light was found as the liquid 
passed into the supercooled region when viewed between crossed nicols. 

Apart from a few determinations by White and Twining (1913) of the viscosity 
of supercooled water (an unsuitable substance as its temperature variation of 
viscosity differs markedly in form from that of the vast majority of liquids),.the only 
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systematic determination on a supercooled liquid has been made on liquid white 
phosphorus by Dobinski (1934). Although the accuracy of these determinations. 
is not great (about 2%), he obtains a curve similar to Figure 3, but his results. 
indicate that the departure from the linear law occurs some 5° c. above the melting 
point, which conflicts with our results. When it is borne in mind that phosphorus. 
can exist in various allotropic forms and that on supercooling the clear liquid 
phosphorus becomes turgid and opaque (a phenomenon which is absent in phenyl 
ether), it is clear that the mechanism of the viscosity anomaly in phosphorus is. 
different from that for phenyl ether discussed here. 

The viscosity measurements are being extended to other liquids, including 
liquid gallium, which readily supercools-and whose structure is comparatively 
simple compared with diphenyl ether, which was chosen initially purely because of 
its convenient melting point. Curves for logy against 1/T are plotted in Figure 5 
for m-chlor-nitrobenzene (melting point 44-2° c.), the point of intersection of the: 
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Figure 5. 


two straight lines being at 44-28° c. The deviations of all points from the straight 
line through points taken at temperatures above the melting point are plotted in 
the insert to Figure 5. It is hoped also to carry out determinations of the surface 
tension and of the refractive index into the supercooled region. 
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DISCUSSION 


Dr. R. Firtu. I think that the most interesting results of Dodd and Hu can be 
-explained on the basis of the “‘ hole ”’ theory of the liquid state (Firth et al. 1939, Firth 1940, 
1941). According to that theory the structure of a liquid is characterized by the presence of 
holes whose size distribution depends on temperature and pressure. The bulk of the sub- 
‘stance outside the holes is supposed to have practically the same structure as the corre- 
sponding solid. Viscosity is considered to be the outcome of the transfer of momentum by 
the Brownian movement of the holes, and its magnitude is determined by the mean free path 
-of the holes, which in turn depends on their mean lifetime. It can be shown that as a result 
the viscosity coefficient 7 contains the Andrade factor exp ¢/RT, where ¢ is the molecular 
‘work function for the ‘‘ evaporation” of a molecule into a hole. The total energy 4U for 
the formation of all the holes in a mol of the liquid is equal to¢Z, where Z is the total number 
‘of molecules at the hole surfaces. 

In a supercooled liquid the distribution function of the holes is essentially different from 
that above the normal melting point, indicating an excess number of very small holes. 
Accordingly Z must drop suddenly when the melting point is passed. On the other hand, 
as no latent heat is released, 4U must remain constant, from which it follows that ¢ must be 
larger for the supercooled liquid than for the normal one. This is exactly what the experi- 
ments of Dodd and Hu have shown. Furthermore, it follows from the theory that the 
‘decrease in Z should be of the order of magnitude of 10%. Hence the increase in ¢ should 
be also of the same order, which agrees with the experimental facts. 
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DELLE RS lO RikEgee DITOR 


The Radiation Pattern of Retinal Receptors 


I have recently made use of microwave theory in order to explain the directional 
“sensitivity of the retina (Toraldo di Francia 1948, 1949). As regards the propagation of 
light in their interior, retinal receptors must be thought of as dielectric waveguides, while 
the more general problem of reception must be dealt with by means of the theory of antennas. 
Each elementary receptor is a microantenna, and its radiation pattern, which, on account of 

the reciprocity theorem (Schelkunoff 1943, p. 478), is the same in emission as in reception, is 
‘responsible for the directivity known as the Stiles and Crawford effect. I pointed out that 
‘the absence of the effect for the rods could be accounted for by their diameter being smaller 
‘and their radiation pattern wider than those of the cones. If the principal lobe of the rod 
‘pattern had a larger angular aperture than the pupil, the directivity of the rods could have 
remained concealed. 

I am indebted to Dr. W. S. Stiles for drawing my attention to a recent paper (Flamant 
and Stiles 1948) where some evidence is given regarding the existence of a directivity of the 
‘rods. Making use of the fact that for a test stimulus of short wavelength, which is viewed 
parafoveally, the absolute threshold of the cones greatly exceeds that of the rods, these 
authors were able to separate the cone diagram from the rod diagram with certainty. The 
results obtained by these two observers (F.F. and W.S.S.) are shown in Figure 1, where the 
upper curves relate to the rods and the lower ones to the cones. It appears that the rods 
‘have almost no directivity. 
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Figure 1. Directional sensitivity curves for two normal Figure 2. Directional sensitivity 
observers. Upper curve: rods. Lower curve: curves for an observer hay- 
cones. (From Flamant and Stiles 1948, by kind ing the maximum sensitivity 
permission of the authors.) point displaced from the pupil 


centre. (From Flamant and 
Stiles 1948, by kind per- 
mission of the authors.) 


However, the same authors have published the curves obtained with the observer R.G.B. 
‘(Figure 2), who presents in a high degree the not rare phenomenon (Stiles 1939, Flamant 
1946, Le Grand 1948) of the point of maximum sensitivity being displaced from the pupil 
‘centre. ‘This is a lucky circumstance for the researcher, since it allows the radiation pattern 
to be measured up to a wider angle than usual. In fact the upper curve of Figure 2 gives 
-evidence in favour of the hypothesis that the rods too have a directivity which in ordinary 
-cases is masked by the remarkable angular width of the major lobe. 
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It seems to me that Figure 2 shows another extremely interesting fact. Almost every 
antenna pattern of high directivity presents, besides the major lobe, one or more side lobes of" 
minor importance, and I had been wondering for some months whether the side lobes of the 
cones could be revealed by a suitable experiment. The answer to this question is given 
again by Figure 2. From the data of the lower curve I have been able to plot the sensitivity 
of the cones against the angle of incidence with respect to their axis. The polar diagram. 
thus obtained is shown in Figure 3, where the curve has been completed by symmetry on the- 


Figure 3. Polar radiation pattern of the cones, Figure 4. Polar radiation pattern of a dielectric- 
deduced from Figure 2. antenna radiating from one of its ends. 


nasal side. A part of the first side lobe is clearly visible. It is rather striking to compare 
this radiation pattern with that of a dielectric antenna, shown in Figure 4. 

It would seem advisable in the future to give all the directional sensitivity diagrams of 
retinal receptors in polar coordinates, as is done in Figure 3 and in all standard textbooks on 
antenna theory. 


Institute of Physics of the University, G. TORALDO DI FRANCIA. 
Firenze, Italy. 
16th April 1949. 


Friamant, F., 1946, Comm. Lab. Inst. Opt., 2, 47. 
FLamant, F., and Sties, W. S., 1948, 7. Physiol., 107, 187. 
Ler Granp, Y., 1948, Rev. d’Opt., 27, 759. 
ScHELKUNOFF, S. A., 1943, Electromagnetic Waves (London and New York : McGraw-Hill). 
eu eat W. = 1939, ey Roy. Soc. A, 127, 64. 
ORALDO DI Francia, G., 1948, Nuovo Cimento, 5, 589; 1949, F. Opt. Soc. Amer., D4 
Atti Fond. G. Ronchi, 3, 137. ae sgn! mene 


463: 


REVIEWS OF BOOKS 


Chambers’ Six-figure Mathematical Tables, by L. J. Comrie. Volume 1,. 
Logarithmic Values, pp. xxii+576 ; Volume 2, Natural Values, pp. xxxvi 
+576. (Edinburgh and London: W. & R. Chambers Ltd., 1948.) 


42s. each volume. 


Chambers’ Mathematical Tables is probably the best known set of tables (with the 
possible exception of the Board of Education’s four-figure set), and the present two volumes 
are in some sense the successor to it, for it appears that in 1944, the centenary of their 
publication, Dr. Comrie was asked to report on that work with a view to bringing it 
thoroughly abreast of modern requirements, and that the new work is the outcome. 

The bulk of the old book was a 200-page table of seven-figure logarithms and just over 
130 pages of natural and logarithmic trigonometric functions; the remainder of its 454 
pages were occupied with a miscellany of geodetic and navigational tables. The new work, 
even if a successor to the old, is in no sense a revision of it. Dr. Comrie has considered 
first the requirements of the present-day physicist and statistician and provided the tables 
he is likely to need, to an accuracy which seems adequate for most work. The logarithms 
are now to six figures, and occupy 180 pages, but data for obtaining eight-figure logarithms 
are also provided (the old book gave ten figures for argument 1 to 1:06 at intervals of 
0-0005). 

In volume 2, the more interesting of the two, we have, for x=1 to, 1000, the square,. 
cube and fourth power (fifth power up to 100), the reciprocal, the factorial and the prime 
factors as well as logx! and-x?, (10x)?, xt, (10x)* and (100x)? and some powers of the 
reciprocals. ‘The factorials above 200 are here published for the first time. From 1,000 
to 3,400 we have the square, cube and prime factors. There is also a table of the first 
1,540 primes, the largest of which is 12,919; evidently this table is in no sense a com-. 
petitor with the specialist tables for the use of those interested in properties of numbers. 
There are three sets of tables of the circular functions, with arguments respectively in 
degrees and minutes, in radians, and in degrees and decimals of a degree. In each case 
the cot and cosec for small angles are given at a reduced interval, and the usual auxiliary 
functions are provided for the functions where they become excessively large. There is a 
table for conversion from rectangular to polar coordinates. 

The hyperbolic functions are well tabulated, as well as natural logarithms,. the inverse 
hyperbolic (and circular) functions and the Gudermannian with its inverse. 

Beyond these elementary functions, a selection of transcendents has been made, the 
choice falling on the gamma function and the probability integral in its two main forms.. 
Here the range has been most carefully considered so that the danger of a user wanting 
values beyond the arguments tabulated is very small. Where so much is provided it is 
ungenerous to complain, but I think many a physicist would have welcomed some tabulation 
of elliptic functions, if only in the form of the elliptic integrals EF and F. 

This very short account of the volume, which does not do justice to the care taken in. 
the arrangement of the tables to make interpolation easy and to ensure the maximum 
convenience in use, also fails to mention another very valuable feature—the little treatise 
on numerical differentiation and differential equations, which comes near the end. ‘There 
is also a great deal of really valuable information in the introduction, and a very helpful. 
bibliography of larger tables at the end. 

In volume two, the 6- and 8-figure logarithms have been mentioned. There is also a 
table of antilogarithms, a function for which tables to more than four figures have been 
remarkably scarce up to the present. Naturally, too, there are tables of the logarithms 
of trigonometric functions, which, as in volume 1, are given for angles expressed in degrees 
and minutes, or in degrees and decimals, or in radians. The S and T functions are also 
given for all these arguments and for seconds of time as well. There are tables of log I(x); 
and of the logarithms of the hyperbolic functions. 
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Taking the two volumes together, it seems reasonable to suppose that they will suffice 
for the everyday needs of most workers who have calculations to do, though each will supple- 


-ment them with tables for his own special work; perhaps he will want a selection of tables 


for statisticians, or tables of Bessel functions, or maybe of factors. From this point of 


view the collection is not the successor to Chambers’ Tables, but to such collections as 


f 


those of Hoiiel or Potin, Barlow, Hutton, Dale or Dwight. Certainly, with Barlow’s tables 
on one side of the new book and Jahnke and Emde on the other, any worker would have a 


-fine collection of tables at hand. 


From a footnote to the preface it appears that there is also an abridged edition of 
these 6-figure tables, in one volume, but the exact contents are not known to the reviewer. 


There is also a suggestion that the 64-page book of four-figure tables prepared by the same 


compiler for Messrs. Chambers is a miniature version of the present two volumes. It 
is certainly a very handy volume, but to the present writer it would have seemed that a 
nearer analogue is Standard Mathematical Tables by Milne-Thomson and Comrie. 

1b: ae 


Some Aspects of Oculo-Refractive Technique, by MALCOLM CHOLERTON. Pp. 148. 


(London: Hammond, Hammond & Co. Ltd., 1948.) 21s. net. 


Since numerous physical, physiological and psychological factors are involved in even 


‘the simplest visual act in any visual experiment or test it is extremely difficult for the 
-observer to rid himself of the psychological elements of experience and memory and to 


.Snow is always “‘ white ’ 


extricate the simple sensation from the complex mental precept which he experiences. 
> and coal is “‘ black’ whatever may be the intensity and spectral 


-composition of the light by which they may be illuminated. When experimenting on vision 
\we are dealing with the consequences of the interaction between a physical agency, radiation 


and living organisms ; and running parallel with these physiological actions, consciousness 
is involved. The analytical methods of physical science, in which the subject of enquiry 


-can be resolved into elementary concepts such as atoms and electrons, cannot be applied 


or can be applied only in a modified form. 
Because of the confusion of all these factors, progress in visual science has been much 


less impressive than in physical science ; and its literature is frequently misleading. And 


how difficult it is to define a “‘ normal observer”’ and to know whether such a one is operating ! 


_Not only may an observer classified as normal be temporarily upset physiologically, but he 


may, perhaps for a reason over which he has no control, be incapable of the normal interest 


-In or attention to the work in hand, or of eliminating some psychological factor that is 


vitiating the results. 

These difficulties may be present during clinical testing by the optical practitioner. 
Fortunately, however, such is the adaptability of the eye and visual mechanism, they do 
not obtrude unduly in the majority of cases. The physical tests which the practitioner 
applies, with such instruments as the ophthalmoscope to search for possible signs of disease, 
and with lenses and physical appliances to measure and treat ametropia and oculo-motor 
imbalance, usually provide acceptable comfort. 

To meet those cases in which the orthodox clinical methods prove to be inadequate, 
various methods are practised by enterprising practitioners seeking to improve their service 


to the community, but they are empirical methods and often unsuccessful. More rational 
; methods are not yet possible because our knowledge of the visual process and of the nature 
- of vision, particularly binocular vision, is too scanty. 


There is much discussion of practical clinical problems amongst opticians, but in the 


-absence of the needed basic knowledge of visual science such discussion amounts to little 


more than dialectics in which the terms used are not precisely defined ; and little or no 


‘progress is made. Investigation is sorely needed by those who are skilled in the methods 


- of research and who will pay heed to the biological aspects of the problems. 


The book under review is concerned with certain clinical methods which will be of 
interest to established optical practitioners. It does not attempt a really critical review 


of such methods but sets out, to quote from the author’s preface, “ to bring together, and 
discuss generally, the various phases of the clinical examination of the eye and visual 


—— 
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function in so far as they extend and should concern the optician”’. Chapters are devoted 
to general routine, anisometropia and aniseikonia, presbyopia, dynamic retinoscopy, 
external and internal examination, visual fields, colour vision testing, tinted lenses, orthoptics, 
secondary correspondences and lighting practice. 

Under several of these headings, however, it is not discussion that is needed, but experi- 
ments and yet more experiments. The author of this book draws attention to the prevailing 
confusion of terminology, particularly concerning dynamic retinoscopy, but unfortunately 
his own style of presentation is, to this reviewer at least, frequently obscure. 

For example, in the preface the author expresses his dissatisfaction ‘‘ with the efficacy 
of our present clinical approach to the analysis and rehabilitation of those forms of visual 
affliction or implication not wholly amenable to solution by the present means at our disposal”. - 
He claims that opticians “are skilled in the physics of the situation, in the differential 
weighting of ametropia and anomalies of binocular vision in terms of the relative symptom- 
potential factor..... ”. In dealing with the fusional reserves in the chapter on general 
routine, and having referred to the blur or break point and negative vergence reflex phase, - 
he says, ‘‘this involves consideration of the aspect whether a blur point can be demonstrated 
by the survival of the negative vergence reflex to a stage which can impose upon the accom- 
modation a preceptual disjugacy with the distance of the target, by extending it reflexly 
further in association with the still surviving negative vergence reflex, until the latter is 
finally overwhelmed.” 

Although this style, which is general throughout the book, is difficult to follow, practising 
opticians will find much to stimulate them. H. H.E. 


A Textbook of Radar, by the Staff of the Radiophysics Laboratory, C.S.LR., 
Australia. Pp. vi+579. (London: Chapman and Hall, Ltd., 1948.) 50s.. 


The relaxation of military secrecy since the end of the last war has brought forth a 
spate of books on radar and related subjects. ‘The present book is the work of the wartime 
radiophysics team in Australia. It was first published in Australia and it is now being made 
available to British readers through Messrs. Chapman and Hall, Ltd. Like the earlier 
American and Italian books, it shows evidence of having been written in the latter half of 
1945. For instance no references are given to the important papers on radar published in 
the Journal of the Institution of Electrical Engineers during 1946 or the corresponding American 
publications. This is a pity for a book published in 1948, and it will limit the value of the 
book for the serious student of radio. 

There are twenty chapters in the book and twenty-one contributors. The editing has 
been carefully done, however, and in consequence there are few obscurities and few definite 
mistakes. The first two chapters deal with the fundamentals of the subject and cover the 
usual subjects—special features of radar, components of a radar set, the echoing process, 
the free space radar equation, normal propagation, super refraction, atmospheric absorption, 
and atmospheric scattering. "The explanations are clear and concise and should be easily 
understood by anyone with a pre-war knowledge of radio. 

Chapters III, IV and V deal respectively with the magnetron, triode power oscillators, 
and modulators. Both the theoretical and the practical aspects of these devices are discussed 
and examples of design problems are included. 

The following four chapters continue the story on radar components and deal with 
transmission lines, waveguides, aerials and aerial duplexing. They occupy 210 pages, and a 
very full discussion of the basic principles and the design of these components is given. 
The examples mentioned cover a wide field. ‘Thus, in the chapter on aerials, there are 
sections on lens aerials, metal plate lenses, electromagnetic horns, Foster scanners and Musa 
systems as well as the more orthodox types. There are, however, one or two omissions. 
Thus little is said about the important subject of back-to-front ratio, and as the ‘aerial 
chapter includes a detailed discussion on the theory of linear arrays it would have been 
helpful to have corrected the widely believed impression that the power gain of aerial 
system is theoretically limited by its size. However, in spite of these omissions these four 
chapters will be of real value to the student. 
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The later chapters in the book (X—XV) deal respectively with the characteristics of the 
radar receiving system—the receiver, the timing unit and the presentation system. Separate 
. chapters are included on local oscillators, frequency converters and amplifiers, and there are 
chapters on display circuits and automatic ranging circuits. For the most part the authors 
have concentrated on giving a dissertation on the principles of circuit design as applied to 
these components, but several examples of complete circuits are given with full explanations. 
The final section of the book (Chapters XVI—XX) discusses radar systems. ‘This is 
prefaced by an introductory chapter dealing with functional and operating characteristics of 
radar systems. Bearing in mind the length of the book, this chapter could have been rather 
‘longer with advantage. The later chapters deal with ground radar, shipborne radar, air- 
borne radar and radar navigation respectively. All the main types of equipment are men- 
-tioned and military as well as civil uses are discussed. ‘The treatment is concise and for 
-the most part accurate. The authors have made one or two small factual errors in describing 
British radar equipment, but these are only of historical significance. An amusing mistake 
-occurs on page 515, where the authors, in describing the CH system, state that the centres 
-of both the upper and lower transmitting and receiving aerials were at the same height. 
Actually they weren’t, but they should have been ! 
The book is well printed and bound and should be of value to any serious student of 
rradio. DENIS TAYLOR. 
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ABSTRACTS FOR SECTIONS 


A Recording Spectrometer for Raman Spectroscopy, by C. H. MILter, D. A. Lone, 
L. A. Woopwarp and H. W. THompson. 


ABSTRACT. A description is given of a photoelectric recording spectrometer for Raman 
spectroscopy incorporating a high-intensity excitation system and a novel amplification 
device. Operation of the exciting mercury arc lamps from 50 c/s. mains imparts to the light 
“signal, and consequently to the current from the photo-multiplier tube, a 100 c/s. periodicity. 
After amplification by a narrow band-pass amplifier, rectification is achieved in a push-pull 
phase-sensitive (homodyne) rectifier, by “‘ mixing ”’ the signal with a fixed voltage of identical 
periodicity. The rectified signal is fed through a filter to a cathode follower valve which 
‘provides the power to operate a recording milliammeter. This method of rectification, in 
conjunction with appropriate filter, effectively reduces the bandwidth of the amplifier, with 
consequent great improvement in the signal/noise ratio. With this arrangement cooling 
‘of the photo-multiplier is unnecessary for normal work. Using a spectrometer of small 
aperture and resolution and a Raman tube of 35 ml. capacity, traces of the Raman spectrum 
-of carbon tetrachloride show no appreciable background noise. 


The Neutrons Emitted in the Disintegration of Lithium by Deuterons, by L. L. 
GREEN and W. M. Gisson. 


ABSTRACT. "The photographic plate method has been used in an investigation of the 
neutrons from the reaction 
*Li+?D > {Be+in. 


The results are in general agreement with those of Richards, but show the structure of 
the neutron groups more clearly. In particular, the shape of the broad group attributed 
‘by Richards to a level in *Be at 3:0+-0°4 Mev. is discussed in detail. 

The total energy release is 15-0-+0-15 Mev., and evidence is put forward for the existence 
of energy levels in *Be at 2:8 Mev., 4:05 Mev., 4:9 Mev., and 7:5 Mev. 
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The Basis of the Electron Theory of Metals, with Special Reference to the Transition 
Metals, by N. F. Morr. 


ABSTRACT. It is shown that the collective electron and London-Heitler models are 
not to be regarded as different approximations to the same exact wave function for solids in 
which, according to the former model, there is a partially filled zone of energy levels. It 
can thus be shown why nickel oxide in the pure state is a non-conductor, although it 


contains an incomplete zone. The properties of the metals nickel, palladium and platinum, 


are discussed in the light of these results ; platinum differs from nickel in that the orbital _ 
contribution to the moment of the elementary magnets is not quenched. A discussion is. — 


given of x-ray absorption edges, and it is shown why exciton lines are absent for metals.. 


The Experimental Basis of Electromagnetism—Part III: The Magnetic Field, by 
N. R. CAMPBELL and L. HARTSHORN. 


ABSTRACT. The principles outlined in previous parts, published in the Proceedings of 


the Physical Society in 1946 and 1948 respectively, which dealt with the direct current circuit 
and with electrostatics, are here applied to magnetism, with the object of showing how the 
basic concepts are defined in terms of the operations actually performed in measuring them. 
This part is confined to a discussion of the vector B, which is shown to be measurable 
independently of a knowledge of any other magnitude. It is solenoidal and can therefore 
be determined at all points, even those within solid bodies. 


The Experimental Basis of Electromagnetism—Part IV: Magnetic Materials, by — 


L. HARrTSHORN. 


ABSTRACT. This paper concludes the inquiry which was undertaken by N. R. Campbell — 


and the author with the object of elucidating the connection between the concepts and. 
principles of electromagnetism and the experimental operations actually performed in the 
Laboratory, that constitute their factual basis. In Part III the vector B was established as 
measureable everywhere, even within solid bodies. The vector H and the scalar .=B/H 
are now shown to be measurable in special circumstances by means of the magnetometer 
and permeameter, but in general their values depend on a hypothesis, which is stated. The 
significant facts concerning the magnetic properties of real materials are briefly reviewed,,. 
and the two-fold aspect of the science is emphasized: the self-consistent mathematical 
theories based on postulates, on the one hand, and the complex of experimental laws on the 
other. It is important that the relation between the two shall be clearly established, and not 
merely assumed as self-evident. . 


Temperature Variation of the Electrical Properties of Nickel Oxide, by R. W. 
WRIGHT and J. P. ANDREWSs. 


ABSTRACT. An experimental method of preparing specimens of NiO in the form of 
coherent strips is described. Simultaneous measurements of electrical conductivity o and 
Hall effect R are then made at various temperatures up to about 700° c. and the thermo- 
electric power dE/dT measured immediately afterwards, all on one sample. The experi- 
ments are then repeated on different samples, and they show that log o, log R, and dE/dT 
are linear functions of 1/7’, where T' is the mean absolute temperature of the specimen. The 
graphs for log o and log| R| against 1/T each consist of two rectilinear portions, one covering 
the high, the other the low temperature results. ‘Thermo-electric measurements yielded 
results only for the high temperature part. 

The theory of impurity semiconductors is developed and found adequate to explain the 
variation, at any rate at higher temperatures. It is found that an activation energy of about 
2 ev. at high temperatures, and from 0-3 to 0-64 ev. at low temperatures, is required for 


excitation of electrons to the conduction band. The concentration of impurity centres’ 


operative in the high temperature range is about 8-6 x 102°, the electrons having a mean free 
path of about 10-*cm. ‘The concentration of free “ carriers’ varies from 10" to 1018 as 
the temperature rises. NiO behaves like a defect semiconductor, ‘‘ hole’? conduction 
predominating. ‘The significance of these results is discussed. 
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